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This  report  summarizes  the  research  performed  during  the  period  1  June  1984 
to  31  May  1985  on  adaptive  control  techniques  for  Large  Space  Structures  (LSS). 
The  research  effort  concentrated  on  two  areas:  (1)  a  further  development  of  an 
earlier  idea  on  how  to  perform  on-line  robust  design  from  identified  models  - 
what  is  referred  to  here  as  adaptive  calibration;  and  (2j  an  analysis  of  slow- 
adaptation  for  adaptive  control  of  LSS.  The  reoort  summarizes  the  results 
obtained  in  these  areas  and  also  includes  Appendices  which  contain  detailed 
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SECTION  1 

BACKGROUND  AND  PROBLEM  STATEMENT 


1.1  HISTORY  OF  CURRENT  RESEARCH  PROGRAM 

'~^The  Large  Space  Structure  (LSS)  research  program  was  originally  for¬ 
mulated  in  late  1982  in  response  to  the  increasing  concern  that 
performance  robustness  of  Air  Force  LSS  type  systems  would  be  inadequate 
to  meet  mission  objectives.  In  particular,  uncertainties  in  both  system 
dynamics  and  disturbance  spectra  characterizations  (both  time  varying  and 
stochastic  uncertainty)  significantly  limit  the  performance  attainable 
with  fixed  gain,  fixed  architecture  controls.  Therefore,  the  use  of  an 
adaptive  system,  where  disturbances  and/or  plant  models  are  identified 
prior  to  or  during  control,  gives  systems  designers  more  options  for  mini¬ 
mizing  the  risk  in  achieving  performance  objectives. 

The  aim  of  adaptive  control  is  to  implement  in  real-time  and  on-line 
as  many  as  possible  of  the  design  functioi^^ow  performed  off-line  by  the 
control  engineer;  to  give  the  controller  ’^telligence"^^  To  realize  this 
aim,  both  a  theory  of  stability  and  performance  of  such  Inherently  non¬ 
linear  controls  is  essential  as  well  as  a  technology  capable  of  achieving 
the  implementation.  As  has  been  noted  by  Astrom  ( 1 981^ recent  advances  in 
each  of  these  fronts  have  brought  us  to  the  position  where  adaptive  con¬ 
trol  has  been  applied  to  many  processes  and  standard  adaptive  controllers 
can  be  purchased,  e.g.,  Egardt  (1981)).  _ 


First  Year  Technical  Results  (1  June  1983  “  31  May  198))) 


Motivated  by  the  uncertain  order  of  LSS  models,  the  first  year  re¬ 
search  was  originally  directed  toward  using  variable  order  adaptive 
systems.  Early  in  the  research,  however,  lack  of  a  well-developed  robust¬ 
ness  theory  for  adaptive  mechanizations  required  a  reexamination  of  the 
problem  at  a  more  fundamental  level,  i.e.,  development  of  model  and  dis- 


turbaiice  uncertainty  bounds  for  which  adaptive  algorithms  would  exhibit 
(stable)  desired  performance.  Adaptive  theory  at  that  time  only  con¬ 
sidered  model  and  disturbance  parameter  uncertainty,  and  hence,  did  not 
provide  the  means  to  obtain  the  relation  between  achievable  performance 
and  unmodeled  effects,  e.g.,  unmodeled  residual  modes  and  small  broadband 
disturbance.  The  theory  required: 

(1)  a  complete  parametric  model  (no  unmodeled  dynamics/disturbance) 

(2)  known  model  order 

(3)  known  high  frequency  gain 

Under  these  conditions,  it  was  possible  to  show  that  the  adaptive  system 
was  globally  stable,  i.e.,  all  the  signals  remained  bounded  for  any  init¬ 
ial  set  of  parameter  values  and  any  bounded  reference  command.  If,  in  ■ 
addition,  the  reference  input  was  "sufficiently  rich",  then  the  parameters 
converged  exponentially  fast  to  a  unique  tuned  set  of  values  (the  paper  by 
Goodwin,  Ramadge,  and  Caines  (1980)  summarizes  the  above  theory). 

For  the  LSS  problem,  it  is  not  possible  to  satisfy  a  single  one  of 
the  above  theoretical  requirements.  In  the  first  place,  the  LSS  is 
theoretically  of  infinite  order  and  a  complete  parametric  model  would  be  a 
partial  differential  equation,  not  the  ordinary  differential  equation  re¬ 
quired  by  the  theory.  Secondly,  and  most  importantly,  the  theory  could 
not  tolerate  even  small  deviations  from  the  above  assumptions.  It  was 
shown  in  simulations  by  Rohrs  et  al .  (1981,  1982)  that  unmodeled  high  fre¬ 
quency  dynamics,  together  with  a  small  high  frequency  noise  -  both  well 
outside  the  controller  bandwidth  -  could  cause  a  parameter  drift  which 
eventually  destabilized  the  system.  This  drift  phenomena  was  also  re¬ 
ported  to  occur  in  on-site  process  control  applications  where  it  was 
necessary  to  re-set  the  parameters  every  so  often  (Wittenmark  and  Astrom, 
1982).  Simulation  studies  involving  LSS  exhibited  similar  problems  where 
the  parameters  would  not  converge  unless  the  initial  parameter  values  were 
close  to  the  unknown  tuned  values  (Gupta,  Lyons  et  al . ,  1981;  Sundararajan 
and  Montgomery,  1982). 


Because  of  these  difficulties,  the  first  year  effort  focused  on  re¬ 
laxing  the  requirements  on  adaptive  theory  by  emphasizing  local  rather 
than  global  results.  The  distinction  between  global  and  local  is  that  in 
local  theory  there  are  restrictions  on  the  magnitudes  and  frequency  con¬ 
tent  of  the  Inputs  as  well  as  the  initial  parameter  values.  In  global 
theory  there  are  no  such  restrictions.  Thus,  the  local  theory  is  more 
practical  and  is  able  to  use  a  priori  information  that  is  available. 

To  summarize,  there  were  two  major  accomplishments  in  the  first  year 
effort; 

(1)  Development  of  a  local  theory  of  adaptive  control  with  broad 
applicability.  These  results  are  reported  in  several  joint 
papers  by  Kosut  and  co-workers  from  1982-198^4  (see  References). 

(2)  Methodology  development  for  LSS  based  on  the  above  local 
theory.  This  result  used  key  ideas  from  parameter  estimation 
and  robust  control  design  under  "slow"  adaptation. 

By  "slow"  we  mean  that  there  is  sufficient  time  to  run  batch  identifica¬ 
tion  before  the  control  system  is  modified.  The  use  of  slow  adaptation  is 
anticipated  for  a  large  class  of  LSS  missions  which  have  quiescient 
periods  useful  for  "calibration".  The  methodology  we  have  developed 
provides  a  guaranteed  level  of  performance  given  an  "identified"  model  of 
the  system  together  with  the  model  error  between  the  system  and  the  iden¬ 
tified  model.  In  fact,  our  methodology  generates  performance  vs.  model 
error  tables  (to  be  stored  in  the  computer)  from  which  the  control  design 
is  determined  strictly  on  the  basis  of  model  error  and  performance  demand, 
rather  than  trial  and  error.  These  results  are  summarized  in  our  annual 
report  (ISI  Report  443,  1  9844)  and  in  Kosut  and  Lyons  (  19844). 

The  status  of  adaptive  theory  at  the  end  of  the  first  year  effort  is 
summarized  in  Table  1-1.  The  maturity  of  adaptive  control  theory  is  com¬ 
pared  to  that  of  linear  control  theory  with  respect  to  the  impact  of 
modeling  assumptions  and  modeling  errors.  Table  1-1  is  primarily 
applicable  to  linear  finite  dimensional  systems.  Adaptive  theory  is  vir¬ 
tually  undeveloped  either  for  »  -  dimensional  linear  systems  or  general 
nonlinear  systems. 
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Second  Year  Results  (1  June  1984  -  31  May  1985) 

The  second  year  effort  continued  with  theoretical  work  on  further 
developing  the  local  theory  with  particular  regard  to  the  fact  that  LSS 
modes  are  uncertain,  densely  packed  and  very  lightly  damped.  These 
characteristics  of  LSS  dynamics  are  among  the  primary  difficulties  faced 
by  practical  identification/ adaptation  algorithms.  A  major  task  in  this 
effort  was  to  investigate  the  properties  of  adaptive  systems  incorporating 
algorithms  with  multi-rate  (two  time  scale)  structures,  and  persistent 
excitation.  These  algorithms  were  studied  under  slow  adaptation.  This 
made  it  easier  to  assure  achievable  specified  performance  levels  despite 
unmodeled  dynamics  and  disturbances.  The  effect  of  speeding  up  the  adap¬ 
tation  has  been  proposed  for  study  in  our  current  work  under  Contract 
FH9620-85-C-009A. 


New  Results  and  Directions 

An  important  series  of  developments  took  place  during  this  second 
year  which  will  have  a  positive  effect  for  our  present  and  future  re¬ 
search.  Our  initial  efforts  at  developing  a  local  stability  theory  have 
just  recently  been  further  advanced  by  Riedle  and  Kokotovic  (19814).  Their 
results  indicate  that  by  combining  the  averaging  theory  of  Hale  (1969) 
with  the  local  stability  theory  of  Kosut  and  Anderson  (1983) i  a  sharp 
stability-instability  boundary  can  be  obtained  for  the  case  of  slow 
adaptation/identification.  Insofar  as  LSS  technology  is  concerned  this  is 
a  major  breakthrough,  because  as  mentioned  before,  slow 
adaptation/identification  is  sufficient  in  many  cases.  This  new  insight 
made  it  very  important  to  more  thoroughly  develop  the  theory  during  this 
second  year.  A  heuristic  explanation  of  our  local  stability  theory  and 
the  method  of  averaging  is  provided  in  Section  2.1.  A  more  detailed  over¬ 
view  is  contained  in  the  summary  paper  by  Kosut  in  Appendix  B. 

Collaborative  Research  Effort 

An  important  point  to  be  made,  regarding  this  research,  is  that 
there  is  a  great  deal  of  collaborative  effort  involved  among  several  re¬ 
searchers  who  share  a  common  interest  in  this  field.  In  fact,  quite  soon 
after  the  initial  work  on  input-output  stability  theory  of  adaptive  sys¬ 
tems  was  published  by  Kosut  and  Friedlander  (1982,  1985)*,  other  pockets 
of  research  groups  began  extending  these  results  in  a  variety  of  direc¬ 
tions,  e.g.,  Ortega,  Praly,  and  Landau  (19814),  Riedle  and  Kokotovic 
(19814),  Kosut,  Johnson,  and  Anderson  (1983).  It  became  apparent  that 
these  researchers  should  get  together  for  informal  meetings.  One  was  ar¬ 
ranged  by  Landau  at  the  University  of  Grenoble  in  July  19814  with  Kosut, 
Ortega,  and  Praly  attending.  A  second  was  arranged  by  Kokotovic  at 
Montana  State  University  in  August  19814,  with  most  of  the  above  research¬ 
ers  attending. 

*  Supported  under  AFOSR  contract  F14920-81 -C-0051  . 
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These  small,  rather  intensely  focused  meetings,  have  proven  to  be  a  major 
catalyst  for  advancing  the  research  efforts  in  this  field.  A  third  meet¬ 
ing,  arranged  by  Kokotovic,  occured  in  Nov.  198^4,  at  the  University  of 
Illinois,  the  purpose  of  which  is  to  organize  efforts  on  the  further 
development  and  use  of  the  averaging  theory  of  Hale  (1969).  A  fourth 
meeting  took  place  in  the  summer  of  1985. 

Also  during  this  period.  Dr.  Kosut  received  AFOSR  approval  to  accept 
an  invitation  from  the  Australian  National  University  to  work  there  as  a 
Visiting  Scholar.  He  spent  a  one  month  period  and  worked  in  close  col¬ 
laboration  with  Professor  Brian  Anderson  and  his  collegues.  The  results 
are  reported  in  Kosut,  Anderson,  and  Mareels  (1985)  and  included  here  as 
Appendix  D. 

The  purpose  of  describing  all  those  meetings  here,  is  to  emphasize 
the  collaborative  effort  that  has  been  involved  in  the  development  of  this 
theory,  and  also  to  indicate  the  degree  of  interest  and  excitement  in  the 
adaptive  control  research  community  about  this  endeavor.  Obviously,  with¬ 
out  the  continuing  support  from  the  various  government  agencies,  e.g,, 
AFOSR,  NSF,  etc.,  none  of  this  would  be  possible,  at  least  not  at  this 
pace.  We  are  also  pleased  to  report  that  another  result  of  these  col¬ 
laborations  is  a  forthcoming  book  by  Kosut,  Anderson,  Kokotovic,  et.  al . , 
on  stability  theory  for  adaptive  systems  (MIT  Press,  Spring  1986). 

Other  Related  Research 

This  program  provides  for  the  development  of  theory  which  will  work 
in  synergy  with  other  related  research  activities  both  at  Integrated 
Systems,  Inc.,  and  at  Stanford  University,  where  Dr.  Kosut  has  a  position 
as  Consulting  Processor  in  the  Department  of  Electrical  Engineering.  At 
Integrated  Systems,  Inc.,  these  related  programs  support  development  of 
hardware  architectures  and  associated  hardware  to  provide  practical  adap¬ 
tive  controls  for  real  DOD/NASA  missions.  One  activity,  with  the  Army 
Munitions  and  Chemical  Command  (AMCCOM),  requires  development  of  program¬ 
mable  board  level  processors  to  implement  adaptive  gun  turret  pointing 
controls.  Control  synthesis  and  simulation  is  followed  directly 


by  programming  of  the  real-time  processors,  thus  sharply  reducing  real¬ 
time  system  development  costs.  Another  activity  is  a  more  general 
treatment  of  the  AMCCOM  problem  where  robust  adaptive  mechanizations  for 
complex  systems  such  as  LSS  will  be  implementable  in  a  real-time  system 
level  hardware  architecture.  This  program  was  funded  under  the  AFOSR/SBIR 
office. 

At  Stanford  University,  Dr.  Kosut  meets  regularly  with  Professor 
G.F.  Franklin,  and  several  students,  all  of  whom  are  working  towards  their 
Ph.D.  in  the  area  of  adaptive  control.  The  students  interests  vary  from 
LSS  to  robotics  to  theory.  Two  of  the  students  have  been  working  on  im¬ 
plementing  algorithms  for  use  on  experimental  flexible  systems  in  the 
Guidance  and  Control  Laboratory  in  the  Department  of  Aeronautics  and 
Astronautics  under  the  direction  of  Professor  R.H.  Cannon,  Their  results, 
both  positive  and  negative,  have  provided  a  strong  impetus  in  the  direc¬ 
tion  of  our  research  program  with  AFOSR. 


1.2  RESEARCH  OBJECTIVES 

The  study  being  reported  on  here,  which  is  part  of  our  ongoing  re¬ 
search,  will  extend  adaptive  theory  and  its  application  to  LSS  problems  in 
several  directions.  These  include  the  following: 


(1)  Theoretical  development  -  The  present  emphasis  is  to  merge 
our  present  local  adaptive  theory  with  the  method  of  averaging 
of  Hale  (1969).  First  attempts  involved  slow  adaptation,  since 
this  covers  many  LSS  situations.  Later  on  we  will  examine  fast 
adaptation.  The  theory  developed  here  provides  for: 

(a)  estimates  of  robustness,  i.e.,  stability  margins  vs.  per¬ 
formance  bounds;  (b)  estimates  of  regions  of  attraction  and 
rates  of  parameter  convergence  to  these  regions.  Later  on,  ex¬ 
tensions  of  the  present  linear  finite  dimensional  adaptive 
theory  will  include  nonlinear  and  infinite  dimensional  plants 
and  controller  structures;  and  (d)  extensions  to  decentralized 
systems . 

(2)  Parameter  adaptive  algorithms  -  Assess  the  behavior  of  dif¬ 
ferent  algorithms.  Including:  gradient,  recursive  least 
squares,  normalized  least  mean  squares,  and  nonlinear  observer 
(e.g..  Extended  Kalman  Filter). 
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(3)  Parametric  models  -  Assess  the  Impact  of  model  choices.  In 

particular,  the  effect  of  explicit  and  implicit  model  choices. 
An  explicit  model,  for  example,  is  a  transfer  function  whose 
coefficients  are  all  unknown.  In  an  implicit  model  transfer 
function,  the  coefficients  would  be  functions  of  some  other 
parameters.  Implicit  models  usually  arise  from  physical  or  ex¬ 
perimental  ground  data,  whereas  explicit  models  are  selected 
for  analytical  convenience. 

{^)  Adaptive  nonlinear  control  -  Although  our  early  effort  is  to 
study  adaptive  linear  control,  there  are  many  LSS  situations 
where  the  control  is  nonlinear. 

Consider,  for  example.  Figure  1-1  which  depicts  two  LSS  control  modes: 

(1)  vibration  suppression,  and  (2)  tracking/ slewing.  In  the  vibration 
suppression  situation  a  "slow"  adaptive  algorithm  (sequential 
identification/adaptation)  may  be  sufficient,  whereas  in  the  tracking/ 
slewing  case  a  "rapid"  algorithm  (parallel  identification/adaptation)  may 
be  necessary.  These  two  possibilities  stretch  adaptive  theory  at  both 
ends,  particularly  with  regard  to  convergence  rate  requirements. 

Note  also  that  the  two  adaptive  controllers  in  Figure  1-1  involve 
not  only  different  convergence  rate  requirements,  but  also  involve  dif¬ 
ferent  controller  structures,  i.e.,  linear  for  vibration  suppression  and 
nonlinear  for  tracking/slewing. 

Clearly  many  combinations  of  models  and  control  structures  are  pos¬ 
sible,  and  at  this  point  it  is  not  possible  to  enumerate  all  that  may  be 
relevant.  However,  by  examining  a  few  sensible  choices,  such  as  those 
compatible  with  typical  LSS  objectives,  the  ensuing  theoretical  develop¬ 
ment  will  indicate  the  necessary  modifications  required  in  each  case. 


Figure  1-1.  Adaptive  Structures 
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1.3  MOTIVATION  AND  LONG-RANGE  GOALS 

The  issues  of  performance  sensitivity,  robustness,  and  achievement 
of  very  high  performance  with  low-order  controllers  can  be  effectively  ad¬ 
dressed  using  adaptive  algorithms.  The  need  to  identify  modal  fre¬ 
quencies,  for  example,  in  high-performance  disturbance  rejection  systems 
has  been  shown  in  ACOSS  (1981)  and  VCOSS  (1982).  The  deployment  of  high- 
performance  optical  or  RF  systems  may  require  on-line  identification  of 
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critical  modal  parameters  before  full  control  authority  can  be  exercised. 
Parameter  sensitivity,  manifested  by  performance  degradation  or  loss  of 
stability  (poor  robustness)  may  be  effectively  reduced  by  adaptive  feed¬ 
back  mechani zations ,  Reducing  the  effects  of  on-board  disturbance  sources 
on  the  system  performance  (disturbance  rejection)  is  particularly  impor¬ 
tant  for  planned  Air  Force  missions.  For  these  cases,  adaptive  control 
mechanizations  are  needed  to  produce  the  three- to-f ive  orders-of-magnitude 
reductions  in  line-of-sight  jitter  required  by  the  mission. 

Research  is  essential  to  identify  the  performance  limitations  of 
adaptive  strategies  for  LSS  control  both  from  theoretical  and  hardware 
mechanization  viewpoints.  The  long  range  goal  of  this  proposed  research 
program  is  to  establish  guidelines  for  selecting  the  appropriate  strategy, 
to  evaluate  performance  improvements  over  fixed-gain  mechanizations,  and 
to  examine  the  architecture  necessary  to  produce  a  practical  hardware 
realization.  The  initial  thrust,  however,  is  to  continue  to  build  a 
strong  theoretical  foundation  without  losing  sight  of  the  practical  imple¬ 
mentation  issues. 

Impact  on  Other  Aerospace  Applications 

The  adaptive  theory  being  developed  for  LSS  control  will  spillover 
to  other  aerospace  applications.  Table  1-2  shows  a  comparison  of  aero¬ 
space  system  adaptive  control  applications.  The  LSS  control  problems 
clearly  involve  aspects  of  the  other  disciplines. 

1.14  REPORT  OUTLINE 

In  the  next  section  we  summarize  the  technical  issues  involved  in 
the  adaptive  control  of  large  space  structures.  Various  detailed  techni¬ 
cal  papers  and  reports  are  included  as  supporting  Appendices. 


TABLE  1-2.  A  COMPARISON  OF  AEROSPACE  SYSTEM 
ADAPTIVE  CONTROL  APPLICATIONS 
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SECTION  2 

PROBLEM  FORMULATION  AND  TECHNICAL  DISCUSSION 


The  development  of  a  design  methodology  for  adaptive  control  of  LSS 
involves  many  different  issues.  In  this  section  we  present  a  selective 
discussion  of  the  theoretical  and  practical  issues  that  seem  most 
relevant.  A  more  in  depth  discussion  of  the  overall  LSS  control  design 
problem  is  presented  in  Appendix  C.  The  discussions  there  cover  various 
types  of  control  design  approaches,  including  both  robust  (non-adapti ve ) 
as  well  as  adaptive.  The  major  points,  however,  are  summarized  in  this 
section. 


2. 1  LSS  CONTROL  PROBLEM  SETTING 

Control  Design  Objectives 

Problems  associated  with  vibration  control  and  accurate  pointing  of 
DEW/LSS  systems  typically  involve  a  combination  of  the  following  control- 
performance  objectives. 

(1)  modal  damping  augmentation  to  enhance  transient  settling  or  im¬ 
prove  quasi-static  vibration  propagation  behavior, 

(2)  stabilization  of  the  attitude  control  system, 

(3)  eigenvector  modification  to  reject  narrow  band  steady-state  dis¬ 
turbances,  and 

(i4)  maneuver  load  management  to  minimize  structural  loads  or  modal 
excitation  (transient  or  steady-state). 

Modeling 

The  basis  for  selecting  a  control  strategy  must  include  an  adequate 
description  of  the  relevant  structural  dynamics  together  with  a  descrip¬ 
tion  of  how  system  performance  is  to  be  measured.  Initially, 


continuum  models  were  suggested  as  the  basis  for  proper  system  design 
since  discretization  of  the  model  could  be  postponed  or  eliminated,  Hughs 
and  Skelton,  1981.  Unfortunately,  practical  spacecraft  configurations  do 
not  present  simple  boundary  conditions  or  simple  shapes,  hence,  p.d.e., 
representations  are  nearly  impossible  to  write.  However,  such  continum 
models  have  provided  useful  insight  into  appropriate  discrete  representa¬ 
tions.  Finite  element  models  can  provide  adequate  fidelity,  at  least  over 
the  frequency  range  needed  for  the  control  design  model,  and  are  supported 
with  sophisticated  software  tools  easily  adapted  to  the  needs  of  control 
design,  ACOSS  (1981). 

Two  Level  Control  Architecture 

The  basic  control  architectures  can  readily  be  combined  into  a  two- 
level  control  system  architecture  consisting  of  a  wide-band,  low-authority 
control  (LAC)  and  a  narrow-band,  high- authority  control  (HAC),  see  Figure 
2-1.  The  LAC  introduces  low  damping  (2-10$)  in  a  wide  range  of  modes  for 
maximum  robustness.  HAC  provides  high  damping  and  mode  shape  adjustment 
in  selected  modes  to  meet  performance  requirement. 

LAC  synthesis  principally  Involves  passivity  methods  and  rate  feed¬ 
back  mechanizations,  usually  with  co-located  actuators  and  sensors,  Aubrun 
(1980),  Iwens,  et.  al  (1980). 

HAC  synthesis.  In  addressing  performance  goals  associated  with 
dynamic  wavefront  and  line-of-sight  error  suppression,  requires  high  modal 
damping  and  mode  shape  changes.  Hence,  the  HAC  is  dependent  on  accurate 
narrow-band  models.  For  such  requirement,  it  is  essential  that  control 
design  techniques  manage  both  dependence  on  model  fidelity  and  system  gain 
in  regions  where  model  fidelity  is  poor.  This  has  generally  been  ac¬ 
complished  using  fixed-gain  robust  control  theory,  e.g.,  Kosut ,  Salzwedel, 
and  Emami  (1983).  With  this  architecture  it  is  likely  that  only  the  HAC 
would  be  tuned  by  an  adaptive  system  since  the  LAC  is  inherently  robust. 


At ( u4t  ors 


bt  rut  I  ura 
livii.tmi  s 


Strnkor 


Figure  2-1.  Features  of  Control  Architecture 


Adaptive  Techniques 

In  general,  uncertainties  in  both  disturbance  spectra  and  system 
dynamical  characteristics  limit  the  performance  obtainable  with  fixed 
gain,  fixed  order  controls,  e.g.,  MAC  system.  The  use  of  an  adaptive  con¬ 
trol  mechanization  where  disturbance  and/or  plant  dynamics  are  identified 
prior  to  or  during  control,  gives  system  designers  more  options  for  mini¬ 
mizing  the  risk  in  achieving  performance  benchmarks. 

In  the  case  of  LSS/DEW  systems,  the  performance  levels  are  extremely 
high.  Hence,  it  is  necessary  that  disturbance  and  plant  models  are  ac¬ 
curately  known.  Since  model  data  obtained  from  ground  testing  is  unlikely 
to  sufficiently  match  the  actual  on-orbit  system,  it  follows  that  on-line 
procedures  are  needed  for  identification  and  control. 
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The  generic  properties  of  closed-loop  system  performance  vs.  struc¬ 
tural  parameter  variations  are  depicted  in  Figure  2-2. 
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Figure  2-2.  Performance  Vs.  Strucutral  Parameter  Variations 

(Generic  Properties) 


Assessment  of  Adaptive  Techniques 
Available  Algorithms 

Many  adaptive  control  and  identification  algorithms  exist  for  lumped 
parameter,  finite-dimensional  linear  systems,  e.g.,  Goodwin  and  Sin 
(1984);  Ljung  and  Soderstrom  (1983).  Most  available  algorithms  can  be 
cast  into  the  form  shown  in  Figure  2-3.  for  example,  a  user  could  select 
from  the  following  catalog  of  model  forms,  control  design  procedures,  and 
parameter  adaptive  mechanisms: 


Model 

Control  Design 

Adaptation 

ARMAX 

Model  Reference 

Gradient 

State-Space 

Self-Tuning 

Recursive  Least  Squares 

Pol e-Placement 

Recursive  Max  Likelihood 
Extended  Kalman  Filter 

Ofcitctw^ 


Figure  2-3.  Adaptive  Control  System 


These  schemes  also  differ  in  terms  of  update  rates.  Typically  the  outer 
control  loop  is  at  a  fast  rate,  whereas  the  parameters  from  identification 
are  updated  more  slowly.  Adaptive  schemes  are  generally  referred  to  as 
recursive  if  the  identification  rate  is  a  fixed  multiple  of  the  controller 
rate.  If  identification  is  used  for  occasional  tuning  or  calibration,  the 
scheme  is  referred  to  here  as  adaptive  calibration. 

Consider,  for  example.  Figure  1-1,  which  depicts  two  LSS/DEW  control 
modes:  (1)  vibration  suppression,  and  (2)  tracking/ slewing.  In  the  vibra¬ 
tion  suppression  situation  a  "slow"  adaptive  algorithm  (sequential 
identification/adaptation)  may  be  sufficient,  whereas  in  the  tracking/ 
slewing  case  a  "rapid"  algorithm  (parallel  identification/adaptation)  may 
be  necessary.  These  two  possibilities  stretch  adaptive  theory  at  both 
ends,  particularly  with  regard  to  convergence  rate  requirements. 
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Note  also  that  the  two  adaptive  controllers  In  Figure  1-1  involve  not 
only  different  convergence  rate  requirements,  but  also  involve  different 
controller  structures,  i.e.,  linear  for  vibration  suppression  and  non¬ 
linear  for  tracking/ slewing. 

Direct  application  of  available  algorithms  to  the  LSS  system  is 
restricted  because  of  a  lack  of  theory  regarding  the  system's  robustness 
to  model  error.  In  addition,  it  is  not  known  how  limitations  imposed  by  a 
decentralized  for  pre-selected  control  architecture  will  effect  achievable 
performance.  Those  two  issues  will  now  be  discussed. 

Robustness  to  Model  Error 

The  use  of  available  identification  and  adaptive  algorithms,  which 
are  based  on  finite-dimensional,  linear  models,  on  LSS  systems  introduces 
a  troublesome  source  of  model  error,  i.e.,  high  frequency  unmodeled  struc¬ 
tural  modes,  of  which  there  are  theoretically  an  infinite  number.  Other 
sources  of  model  error  include  uncertain  actuator/sensor  dynamics  and 
neglected  nonlinearities  in  Joints  and  damping  mechanisms. 

Although  available  theory  can  handle  a  finite  number  of  parameter  er¬ 
rors,  it  cannot  deal  effectively  with  other  types  of  model  error, 
specifically,  the  unmodeled  high  frequency  structural  modes  and  dynamics. 
It  was  shown  in  simulations  by  Rohrs  et  al.  (1981,  1982)  that  unmodeled 
high  frequency  dynamics,  together  with  a  small  high  frequency  noise  -  both 
well  outside  the  controller  bandwidth  could  cause  a  parameter  drift  which 
may  eventually  destabilize  the  system.  This  drift  phenomena  was  also 
reported  to  occur  in  on-site  process  control  applications  where  it  was 
necessary  to  re-set  the  parameters  every  so  often  (Wittenmark  and  Astrora, 
1982).  Simulation  studies  involving  LSS  exhibited  similar  problems  where 
the  parameters  would  not  converge  unless  the  initial  parameter  values 
where  close  to  the  unknown  tuned  values  (Sundararajan  and  Montgomery, 


Robustness  to  model  error  is  more  well  understood  in  the  context  of 


nonadaptive  linear  control  theory,  e.g.,  Doyle  and  Stein  (1961);  Zames  and 
Francis  (1983);  Chen  and  Desoer  (1982).  The  common  theoretical  basis  for 
these  robustness  theories  is  the  input-output  view  of  feedback  systems 
(Zames,  1966;  Desoer  and  V idyasagar ,  1975,  Safonov,  1980).  Another  view¬ 
point  on  robustness,  which  follows  more  along  the  lines  of  Liapunov 
stability  theory,  characterizes  the  solutions  of  perturbed  nonlinear  ordi¬ 
nary  differential  equations,  e.g.,  LaSalle  and  Lefschetz  (1961),  Hale 
(1969). 

A  main  intent  of  our  research  program  has  been  to  merge  the  input- 
output  view  and  the  Liapunov  stability  view  with  adaptive  mechanizations 
to  develop  a  theory  of  robust  adaptive  control.  Some  of  the  groundwork 
has  already  been  accomplished  (Kosut,  et.  al.,  1982-1985)  and  as  such,  we 
are  now  in  a  good  position  to  address  those  issues  as  they  related  to  LSS 
systems. 

Decentralized  Adaptive  Control 

Limitations  on  control  authority  and  the  information  pattern  are  the 
main  features  of  the  decentralized  control  problem.  The  general  structure 
of  such  a  decentralized  control  system  is  illustrated  in  Figure  2-^.  The 
dashed  lines  indicate  a  partial  information  exchange ,  e.g.,  the  local  con¬ 
troller  receives  reference  commands  (or  discretes)  from  a  higher  level 
control  (the  coordinator)  and/or  information  from  other  local  controllers 
in  the  form  of  an  "aggregated"  state. 


Figure  2-U.  Decentralized  LSS  Control.  Dashed  Lines  Indicate  a  Partion 
Exchange  of  Information 


For  LSS/DEW  systems,  decentralization  normally  results  because  a 
natural  separation  is  physically  or  geographically  present  between  func¬ 
tional  components  of  the  system.  For  example,  some  decompositions  result 
from  spatial  differences;  weak  dynamic  interaction  effects  can  be  easily 
identified.  A  decomposition  also  occurs  from  temporal  differences; 
phenomena  occurring  at  different  time-scales,  e.g.,  a  separation  between 
fast  and  slow  modes  or  between  low  frequency  and  high  frequency  effects. 
For  example,  groups  of  the  modes  can  be  separately  controlled  by  separate 
controllers  which  do  not  destabilize  each  other.  Specific  combinations  of 
weak  dynamic  coupling  and  separation  of  slow  and  fast  modes  can  often  be 
identified,  e.g..  Figure  2-5. 

A  number  of  very  useful  results  are  available  for  non-adaptive 
decentralized  systems,  e.g.,  Vidyasagar  (1981),  Siljak  (1978).  Those 
results  show  how  total  system  properties  are  dependent  on  subsystem  and 
Interconnection  properties. 
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Figure  2-5.  Weak  Dynamic  Coupling  with  Slow  and  Fast  Modes 


Much  less  is  known  about  adaptive  decentralized  control,  although 
there  are  some  promising  preliminary  results  available,  e.g.,  loannou 
(198M).  There  is,  however,  a  similarity  in  robustness  theory  and  the 
theory  of  decentralized  control  which  we  hope  to  explore  and  exploit. 

Research  Objectives 


The  dream  of  adaptive  control  researchers  has  always  been  to  develop 
a  "global"  stability  theory,  i.e.,  performance  is  guaranteed  independent 
of  initial  conditions  and  disturbance  spectra.  Such  results  have  been  ob 
tained,  but  as  mentioned,  the  available  theory  requires  that  a  finite 
dimensional  parameterization  exists  such  that  an  exact  matching  condition 
is  obtained  (see  e.g.  Ljung  and  Soderstrom,  1983:  Goodwin  and  Sin,  198')). 
Such  a  stringent  requirement  is  impossible  to  obtain  in  an  LSS  system, 
principally  because  there  are  many  sources  of  model  error  which  defy  an 
exact  parameterization  (see  Section  1). 


2 


Because  of  these  difficulties,  recent  research  efforts  have  focused 
on  relaxing  the  requirements  on  adaptive  theory  by  emphasizing  local 
rather  than  global  results.  The  distinction  between  global  and  local  is 
that  in  local  theory  there  are  restrictions  on  the  magnitudes  and  fre¬ 
quency  content  of  the  inputs  as  well  as  the  initial  parameter  values.  In 
global  theory  there  are  no  such  restrictions.  Thus,  the  local  theory  is 
more  practical  and  provides  robustness  to  model  error  by  utilizing  a 
priori  information  that  is  available.  In  addition,  restrictions  imposed 
by  a  decentralized  architecture  can,  in  principal,  also  be  accounted  by 
the  form  of  the  theory. 

The  basic  ideas  for  local  stability  rest  on  two  fundamental  stability 
theories.  One  is  the  small  gain  theory  of  Zames  (1966),  and  the  other  is 
the  method  of  averaging  as  described  by  Hale  (1969).  These  theories  have 
different  origins  and  to  some  extent  have  application  in  different  regimes 
for  nonlinear  systems.  Small  gain  theory  determines  the  stability  of 
trajectories.  With  regard  to  LSS  problems,  small  gain  theory  is  generally 
applicable  to  "fast"  recursive  adaptive  control;  averaging  theory  is  ap¬ 
plicable  to  "slow"  identification  and  then  control,  i.e.,  adaptive 
calibration.  One  can  envision  small  gain  theory  as  providing  a  macrosco¬ 
pic  view  while  averaging  provides  a  microscopic  view.  Our  intention  is  to 
merge  these  as  much  as  possible  and  thus,  broaden  the  application  of  each 
approach,  particularly  for  the  LSS  system.  This  is  not  without  precedent. 

Small  gain  theory  was  applied  originally  to  continuous-time  adaptive 
systems  by  Kosut  and  Friedlander  (1982,  1985)  and  more  recently  extended 
by  Kosut  and  Johnson  (1984),  and  Kosut  and  Anderson  (1984).  The  method  of 
averaging  was  applied  to  continuous-time  adaptive  systems  with  almost  pe¬ 
riodic  inputs  for  the  first  time  by  Riedle  and  Kokotovic  (1984)  and  Astrom 
(1984).  Extensions  to  more  general  input  classes  have  been  obtained  by 
Kosut,  Anderson,  and  Mareels  (1985).  The  averaging  theory  has  also  been 
able  to  accurately  predict  the  drift  phenomena  observed  by  Rohrs,  et  al . 


Although  small  gain  results  apply  equally  to  continuous  or  discrete¬ 
time  systems,  the  averaging  theory  is  underdeveloped  for  the  discrete¬ 
time  case,  A  major  proposed  early  effort,  therefore,  is  to  develop  a 
discrete-time  averaging  theory  useful  for  sample-data  adaptive  systems. 
This  will  provide  foundation  for  later  efforts  in  sampled-data  adaptive 
systems  with  multi-rate  and  decentralized  processors. 


2.2  SUMMARY  OF  UNDERLYING  THEORY 

In  this  section  we  present  a  heuristic  discussion  which  summarizes 
the  underlying  theory.  As  mentioned,  the  basic  ideas  rest  on  two  fun¬ 
damental  theories:  1)  is  the  small  gain  theory  of  Zames  (1966),  and  2) 
the  method  of  averaging  as  described  by  Hale  (1969). 

To  heuri sti cally  describe  our  approach,  consider  the  feedback  system 
in  Figure  2-6,  below  where  and  G^  represent  operators. 


Figure  2-6.  Feedback  System 


Small  gain  theory  asserts  that  if 


Gain  (G^ )  .  Gain  (G^)  <  1  (1) 

then  the  feedback  system  is  stable.  In  robustness  theory,  the  feedback 
loop  is  arranged  so  that  G^  represents  all  of  the  nominally  zero  uncer¬ 
tain  elements  of  the  system.  The  operator,  G^,  referred  to  as  a  return- 
difference  operator,  is  everything  else.  It  is  always  possible  to  arrange 
a  feedback  loop  this  way  (Safonov,  1978).  Thus,  G^ ,  is  dependent  on  the 
presumed  uncertainty  location.  Hence,  stability  margin  is  the  maximum  al¬ 
lowable  uncertainty  which  guarantees  stability,  i.e.. 


Stability  Margin  =  1/Gain(G^)  (2) 

For  linear-time-invariant  systems  it  is  not  difficult  to  calculate 
Gain  (G^);  in  fact,  this  can  be  accomplished  in  the  frequency-domain. 

These  results  parallel  Bode  analysis,  and  offer  the  engineer  a  very  useful 
design  tool.  At  the  present  time,  however,  no  such  tool  exists  for  adap¬ 
tive  systems.  The  main  difficulty  is  that  in  the  adaptive  case  G^  is 
nonlinear.  Nonetheless,  certain  types  of  nonlinearities  in  G^  can  result 
in  a  frequency- domain  test,  e.g.,  Popov  criterion,  hyperstability  con¬ 
cepts,  positivity,  etc.  These  latter  results  follow  from  passivity 
theorems,  and  interestingly  enough,  this  is  the  main  tool  in  current  adap¬ 
tive  theory  for  proving  stability.  Unfortunately,  however,  by  proceeding 
this  way,  it  can  be  shown  (Kosut,  1982a)  that  the  resulting  stability  mar¬ 
gins  are  very  small  and  are  easily  violated  in  any  realistic  environment. 


The  approach  we  are  proposing  here  is  to  by-pass  the  conservative 
passivity  conditions  and  rearrange  the  adaptive  system  so  that  the  two 
operators  are  intrinsically  small.  This  requires  introducing  the  notion 
of  a  tuned  controller  (Kosut,  1985)  and  a  tuned  set  model  representation 
of  the  plant,  which  accounts  for  uncertainty.  As  a  result  of  this  trans¬ 
formation  we  can  introduce  averaging.  Hence,  a  less  conservative  calcula¬ 
tion  is 

Stability  Margin  =  1/Avg  Gain  (G^  (3) 

The  averaged  gain  is  significantly  smaller  than  the  usual  gain  calcula¬ 
tion,  and  thus,  produces  a  much  larger  stability  margin,  i.e.,  greater 
performance . 

In  Section  2.3  we  outline  this  approach.  Further  details  using  small 
gain  theory  can  be  found  in  Kosut  and  Anderson  (198^)  with  averaging 
results  in  Reidle  and  Kokotovic  (198i»)  and  in  Kosut,  Anderson,  and  Mareels 
(1985).  Appendix  B  contains  a  detailed  report  summarizing  the  applica¬ 
tion  of  averaging  methods  to  adaptive  systems. 

2.3  ADAPTIVE  CONTROL  THEORY 

In  this  section  we  will  discuss  issues  and  new  directions  in  adaptive 
control  theory.  To  illustrate  the  ideas  we  will  consider  a  continuous-time 
generic  representation  for  control  or  identification  of  a  scalar  (single- 
input-single-output)  plant.  In  all  cases  the  comments  apply  equally  to 
discrete-time,  as  well  as  multivariable  systems.  The  generic  adaptive 
system,  shown  in  Figure  2-7,  is  described  by. 


c'  X 

X  =  Ax  +  bz'0 

Dx 

0  =  Yze 

«  p 

where  0(t)  =  0(t)  -  0^  e  R  is  the  parameter  error  vector;  0(t)  is  the 

current  parameter  vector  estimate;  0^^  is  a  constant  tuned  parameter  vector 
P  n 

setting;  z(t)  e  R  is  the  regressor  vector;  x(t)  e  R  is  the  system  state 
consisting  of  plant  and  filter  states;  e(t)  is  the  error  signal;  and  Y  >  0 
is  the  adaptive  gain.  The  signals  e^(t)  and  z, (t),  referred  to  as  the 
tuned  error  and  tuned  regressor,  respectively,  are  outputs  from  the  ideal 
tuned  system  where  0(t)  =  0,  i.e.,  0(t)  =  0^. 
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Figure  2-7.  Generic  Adaptive  System 


Appendix  (A)  shows  how  the  nonlinear  system  (Jj)  arises  in  the  analysis  of 
adaptive  control  or  identification.  More  detail  on  the  development  of 
(2.1)  can  be  found  in  Kosut  and  Friedlander  (1982,  1985),  and  Kosut  and 
Johnson  (19811). 

Global  Stability  Theory 

The  assumptions  required  to  insure  global  stability  and  convergence 
of  (2.1)  are  quite  strict,  e.g.,  Narendra,  Lin,  and  Valavani  (1980),  Good¬ 
win,  Ramadje,  and  Caines  (1980).  These  assumptions  include; 

(A1 )  (t )  — >  0  as  t  — >  “ 

(A2)  z(t)  persistently  exciting  (PE) 

(A3)  c’(sI-A)  ^b  is  strictly  positive  real  (SPR) 

In  all  practical  cases  these  assumptions  are  violated.  Assumption  (A1) 
implies  that  there  are  no  continually  acting  unmodeled  noises  or  distur¬ 
bances.  Assumption  (A2)  can  not  be  guaranteed  unless  (A1)  holds.  The  SPR 
assumption  in  (A3)  is  the  most  restrictive,  since  it  requires  a  priori 
knowledge  of  model  order,  relative  degree,  and  high  frequency  gain,  all  of 
which  are  not  available  or  infeasible  to  obtain.  Moreover,  these  require¬ 
ments  make  no  sense  at  all  for  an  LSS  whose  order  is  theoretically 
infinite.  Unfortunately,  even  small  deviations  from  these  assumptions  can 
critically  disrupt  an  adaptive  system,  e.g.,  Rohrs  et  al .  (1981,  1982), 
Reidle,  and  Kokotovic  (1984). 


Example  of  Adaptive  Calibration 


The  basic  problem  with  control  based  on  identified  modeling  is  that 
without  a  measure  of  model  error  it  is  very  easy  to  destabilize  the  system 
-  particularly  when  the  goal  is  high  performance  -  as  in  an  SBL.  Adaptive 


calibration  is  an  approach  which  incorporates  a  measure  of  model  error 
with  robust  control  design  in  an  iterative  way  so  that  identification  is 
performed  only  where  it  is  needed.  One  proposed  adaptive  calibration  sys¬ 
tem  is  shown  in  Figure  2-8,  with  test  results,  using  the  CSDL  #2  model, 
shown  in  Figure  2-9.  The  adaptive  calibration  procedure  involves  the  fol¬ 
lowing  steps: 


Step  1:  The  model  M(s)  is  a  10-mode  model  which  has  been  obtained 

from  I/O  data. 

Step  2:  Estimate  6(a))  =  model  error  vs.  frequency  using  FFT.  This 

is  dashed  curve  in  Figure  2-9(a). 

Step  3:  Using  the  identified  model  M(s)  and  the  model  error  6,  syn¬ 

thesize  a  robust  control  (e.g..  Appendix  C). 

Step  11-5:  Calculate  “  stability  margin  (or  6  (m)  -  perfor¬ 

mance  margin).  This  is..dark  curve  in  FigBre  2-9(a). 

Compare  to  model  error  6(a)).  Both  plotted  in  Figure 
2-9(a).  If  6(a))<6  (o))  go  to  Step  7  and  implement  control¬ 

ler.  Otherwise,  g8  to  Step  6. 

Step  6:  Modify  filter  windows,  number  of  parameters,  or  input 

spectrum  and  then  repeat  Step  1  to  obtain  new  ID  model. 
Figure  2-9 (b)  shows  result  of  ID  after  one  mode  is  added  in 
frequency  domain  region  where  test  fails. 

Step  7:  Implement  controller. 

These  preliminary  results  are  quite  promising.  A  major  portion  of 
our  initial  research  effort  has  been  to  understand  the  nature  of  such 
schemes  at  a  fundamental  level. 
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Figure  2-8,  A  Procedure  for  Adaptive  Calibration  of  LSS 
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a).  Comparison  of  Allowable  and  Estimated  Model  Errors 


•  PLOTS  SHOW  WHERE  ERROR  IS  >  ALLOWABLE 


Figure  2-9(b).  Data  Before  and  After  ID  Cycle 


Local  Stability  Theory 


Boundedness  Results  (Fast  Adaptation  -  High  Gain) 

With  regard  to  practical  adaptive  design  techniques,  the  need  for  the 
SPR  condition  -  hence,  the  restrictive  modeling  assumptions — can  be 
eliminated  by  considering  conditions  for  local  stability  rather  than 
global  stability,  see  e.g.,  Kosut  and  Anderson  (198^4).  The  term  'local' 
refers  to  the  use  of  known  restrictions  on  the  system  external  inputs,  un¬ 
certain  parameters,  and  unmodeled  dynamics.  For  example,  since  persistent 
excitation  induces  exponential  stability  (Anderson,  1977),  and  since  an 
exponentially  stable  system  is  inherently  robust,  it  is  logical  to  expect 
that  unmodeled  dynamics  could  be  robustly  tolerated.  Other  mechanisms-- 
including  persistent  excitation — can  ensure  stability  of  the  adaptive 
system,  without  SPR,  provided  certain  other  restrictions  are  enforced, 
e.g.,  slowly  varying  signals,  approximate  SPR,  gain  retardation,  and 
restricted  signal  magnitudes  and  bandwidths  (loannou  and  Kokotovic,  1982- 
1984):  Kosut,  1981)).  It  is  our  intention  to  utilize  these  theoretical 
results,  whenever  appropriate,  in  the  adaptive  LSS  study,  as  this  local 
theory  emerges. 


Stability  Results  (Slow  Adaptation-Low  Gain) 

The  above  local  theory,  which  proceeds  from  the  Small  Gain  Theory  of 
Zames  (1966),  and  Desoer  and  Vidyasager  (1975)  can  be  considered  as  giving 
conditions  for  boundedness  (in  the  input-output  sense)  rather  than 
stability  (in  the  Liapunov  sense).  The  recent  work  of  Reidle  and 
Kokotovic  (1981))  utilizes  the  averaging  theory  of  Hale  (1969)  to  obtain  a 
sharp  stability  -  instability  boundary  for  slow  adaptation,  i.e.,  when  Y 
in  (2.1)  is  sufficiently  small.  To  apply  the  averaging  theory  requires 
transforming  (2.1)  to  the  following  form: 


C  =  ae.  +  yg  (t.e.c.Y) 


where  f  and  g  are  time-varying  nonlinear  functions.  As  y  -->  0  the 
stability  of  (2.2)  is  identical  to  the  stability  of  the  "averaged"  system 


e  =  Yf„(e) 

^  =  A^ 


where 


(6) 


fo(Q) 


lim  ^  S  f(t, 0,0,0)  dt 

T  — >  a.  ^  o 


(7) 


Without  going  into  the  details  here,  it  can  be  shown  that  there  exists  a 
sufficiently  small  positive  Y  such  that  (6)  is  stable  —  and  hence,  (5)  is 
stable  -  if  and  only  if  a  persistent-excitation  type  of  condition  is 
satisfied.  The  important  results  here  may  be  summarized  as  follows: 

(1)  For  sufficiently  slow  adaptation,  there  exists  a  sharp  stability 
-  instability  boundary. 

(2)  The  conditions  do  not  require  SPR,  but  do  require  persistent 
excitation. 

Much  work  has  to  be  done  to  further  develop  this  theory  in  general  and 
to  tie  it  together  with  the  boundedness  results  by  Kosut,  et.  al.  For  LSS 
in  particular  it  is  necessary  to  extend  the  theory  to  discrete-time ,  ®  - 
dimensional  and  nonlinear  systems.  Note,  however,  that  the  form  of  the  sys¬ 
tem  (2.2)  required  by  the  averaging  theory  can  be  nonlinear. 


Extensions  to  Discrete-Time 

An  important  step  in  our  research  is  to  develop  the  averaging  theory 
of  Hale  (1969)  to  discrete-time  adaptive  systems.  There  appears  to  be  no 
theory  available  at  the  present  time.  The  obvious  extension  is  to  consider 
the  discrete-time  system  analagous  to  (5),  i.e., 


A9  =  Yf(t,9,C,Y) 

AC  =  AC  +  Yg(t,e,C,Y) 


(8) 


where  A  is  the  finite-difference  operator  (Ax)(t)  =  x(t+1)-x(t)  and  t 
takes  on  discrete-values  t  =  (1,  1,  2,  ...  ).  One  can  make  the  conjecture 
that  the  discrete-time  averaging  theory  states  that  as  Y  — >  Od,  the 
stability  of  (2.2)'  is  identical  to  the  stability  of  the  averaged  system 


A0  =  Yf,(e) 


(9) 


AC  =  AC 


where 


lim  -1  I  f(t,e,0,0) 

T  ->  00  t=0 


(10) 


We  will  conduct  research  to  determine  if  this  conjecture  is  true,  or  what 
the  necessary  modifications  Involve. 


?.U  ADAPTIVE  CONTROL  STRUCTURES 


In  this  section  we  explore  variations  of  (4)  arising  from  different 
models,  control  structure  definitions,  and  parameter  update  laws.  These 
variations  will  be  useful  in  different  LSS  control  applications. 


Structure  of  the  Parameter  Adaptive  Algorithm 


The  adaptive  parameter  update  algorithm  in  (iJ),  0  «  Yze,  although  use¬ 
ful  for  analysis,  is  too  simple  in  practical  cases.  The  more  general  form: 


0  = 


Y  h(t,2,e) 


(n) 


can  arise  from  normalization. 


h(t,z,e)  =  ze/(1  +  |z| 


(12) 


least  squares  considerations. 


h(t,z,e)  =  R(t)  ^  ze/(1  +  z'  R(t)  ^  z) 
R(t)  =  z(t)  z(t)' 


(13) 


dead- zones, 


or  combinations  of  all  of  the  above.  These  all  have  their  use  and  charac¬ 
teristics  properties  [see  Goodwin  and  Sin  (198^))  or  Ljung  and  Sodestrom 
(1983)  for  a  survey].  Which  combination  of  these  or  other  structures  is 
useful  for  LSS  n3eds  to  be  determined, 

Multi-Rate  Adaptive  Control 

Multi-rate  adaptive  schemes  (or  two-times  scale  schemes)  refers  to  an 
adaptive  structure  where  the  adaptive  control  parameters  are  updated  more 
slowly  than  the  basic  control  sampling  rate.  This  scheme  allows  a  period  of 
plant  identification  to  be  followed  by  a  change  in  control  parameters. 
Allowing  more  time  in  the  identification  phase  also  allows  for  more  reliable 
pass-fail  tests.  For  example,  we  can  perform  a  fit-error  test  after  para¬ 
meter  identification  on  the  old  data  as  well  as  estimating  its  information 
content.  If  either  test  fails  then  the  control  gains  are  not  updated.  If 
they  both  pass  then  the  gains  can  be  updated  in  the  direction  of  maximum 
information.  We  can  also  include  a  robustness  test  based  on  a  priori  (or 
new)  model  error  bounds.  Ideas  such  as  these  have  been  considered  for  LSS 
applications  by  Sundararajan  and  Montgomery  (1982), 


Hybrid  Adaptive  Control 

If,  in  the  multi-rate  adaptive  structure,  the  basic  control  sampling 
rate  becomes  continuous,  then  the  structure  is  referred  to  as  a  "hybrid" 
structure:  the  feedback  is  continuous  but  is  updated  at  discrete  times. 


Such  a  scheme  reduces  Intersample  ripple  and  also  tends  to  Increase  band¬ 
width  and  noise  attenuation.  These  structures  may  prove  quite  practical  for 
the  LSS. 


Linear  Finite-Dimensional  Models 

The  modeling  assumptions  influence  the  adaptive  control  structure. 
Assume  that  the  plant  to  be  controlled  can  be  accurately  represented  by  a 
finite-dimensional  model  of  the  form 


B  (u,  y)' 
m  mm  m 


where  u  is  the  input,  y  is  the  output,  a  stable  nxn  matrix,  and 

is  controllable,  (see  e.g.  ch.  2,  Kailath,  1980).  Under  these  conditions, 

the  transfer  function  P  (s,0)  from  u  into  y  is: 

m 


e  ,s"-'  .  e  ,s"‘^  .  . 

P  (S.9)  -  ^ - 

“  s"  .  0  s"''  .....  e 

1  n 


This  type  of  model  is  referred  to  as  an  explicit  parametric  model.  The  ad¬ 
vantage  of  this  representation  is  that  the  parameter  vector  6  appears 
linearly  in  (15).  In  Appendix  A  we  discuss  how  this  form  is  compatible  with 
equation  error  and  output  error  identification  algorithms.  Moreover,  these 
algorithms  can  be  transformed  to  the  form  of  (U)  and  the  theory  discussed 
earlier  can  be  applied.  Also,  in  adaptive  control,  the  linear  relation  be¬ 
tween  y  and  0  in  (15)  simplifies  the  transformation  from  model  parameter 
estimates  0  to  control  parameters,  i.e.,  k  «  f(0). 
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There  are  many  cases  where  the  explicit  model  (15),  although  analyti¬ 
cally  tractable,  is  oblivious  to  useful  information  obtained  from  a  physical 
analysis.  Consider  the  impl ici t  linear  modal, 


y  =  c(a) ' 


X  =  A(a)x  +  b(a)u 
m  m 


where  A(a),  b(a),  c(a)  are  nonlinearly  dependent  on  the  parameter  vector 

l< 

a  c  n,  where  is  the  parameter  space,  a  subset  of  R  .  This  model  (2,9) 
retains  the  physical  meaning  of  the  parameters,  e.g,,  in  an  LSS  the  elements 
with  a  can  consist  of  equivalent  masses,  spring  constants,  and  dampers. 

There  will  be  far  fewer  of  these  parameters  than  there  are  coefficients  in 
the  equivalent  transfer  function,  i.e.,  the  transfer  function  P^(s,a)  from  u 
into  y  for  (17)  is: 


b,(a)3"  '+...+  b„(a) 

- 

s'^  +  a^  (a)s‘’  +  . . .  +  aj.j(a) 


It  is  clear  that  the  parametric  model  of  (18)  can  be  made  identical  to  (15) 
simply  by  equating  each  coefficient  i.e.,  6^  =  a^ (a) ,  6^  =  3^(0),  and  so  on. 
Since  the  model  (16)  is  linear  in  0  whereas  (18)  is  nonlinear  in  a,  there  is 
a  great  difference  in  the  identification  schemes.  For  (18)  there  are  fewer 
parameters  to  be  identified,  but  several  nonlinear  relations.  For  (16), 
there  are  many  more  parameters  to  be  identified  but  all  enter  linearly  into 
the  model.  In  LSS  modeling  it  may  be  prudent  to  examine  the  effect  of  using 
implicit  parametric  models,  because  of  the  reduced  number  of  parameters  to 
be  identified.  This  type  of  model,  however,  leads  to  more  complicated  forms 
than  (iJ)  and  involves  a  nonlinear  observer.  One  such  scheme  using  an  EKF  is 
discussed  next. 


Identification  using  EKF 

The  implicit  parametric  linear  model  (18)  is  compatible  with  an  iden¬ 
tifier  based  on  the  extended  Kalman  filter  (EKF).  The  EKF  parameter 
estimator  has  the  form: 


o  =  (y  -  cCa) ’x) 


X  =  A(&)5c  +  b(a)u  +  l<^(5)(y  -  c(a)'x) 


(19) 


where  the  gain  k^(S)  and  k^(a)  are  the  Kalman  filter  gains  obtained  by 
linearizing  (17)  about  the  current  estimate  x.  Although  EKF  has  been  exten¬ 
sively  studied  (e.g,,  Ljung  and  Soderstrom,  1983)  there  still  remains  no 
concrete  proof  of  guaranteed  convergence  for  identification  and  no  proof  of 
stability  in  the  adaptive  case,  i.e.,  when  a  is  used  to  generate  a  control 
law.  In  as  much  as  explicit  parametric  models  may  better  represent  the  LSS, 
we  will  develop  the  theory  in  this  area.  There  have  been  some  attempts  in 
this  direction  by  Safonov  and  Athans  (1978)  and  Vidyasagar  (1980).  In  this 
latter  work  it  is  assumed  that  certain  signals  are  bounded  in  order  to  prove 
convergence.  In  the  adaptive  case  this  assumption  needs  to  be  proved  first, 
because  the  signals  that  need  to  be  bounded  are  inside  the  adaptive  loop, 
and  hence,  dependent  on  the  parameter  estimator.  The  approach  we  take 

is  to  apply  the  averaging  theory  of  Hale  (1969)  by  transforming  the  EKF 
based  adaptive  system  to  the  form  of  (5). 


Infinite-Dimensional  Systems 


Systems  of  the  form  (11)  not  only  presume  a  linear  finite  dimensional 
plant  but  also  a  linear  finite-dimensional  controller.  Extending  the  theory 
to  the  use  of  infinite  dimensional  plant  and  finite-dimensional  controller 


may  not  prove  to  be  too  difficult,  e.g.,  Kosut  and  Friedlander  (1982). 
However,  the  case  for  a  distributed  parameter  controller  (DPC),  which  is 
perhaps  more  realistic  for  LSS,  is  more  difficult.  What  is  involved  here  is 
adjusting  a  few  parameters  in  the  DPC  from  discrete  spatial-temporal 
measurements.  This  is  certainly  an  area  for  basic  research. 

Nonlinear  Control  Structures 

During  tracking  and  slewing  the  control  can  be  nonlinear,  involving 
switching  curves  (time-optimal)  and  actuator  saturations.  Adaptive  schemes 
here  must  be  rapid  (high  gain)  and  must  account  for  nonlinearities.  Domi¬ 
nant  kinematic  or  structural  nonlinearities  will  also  engender  nonlinear 
control  structures. 

Lattice  Filters  for  Adaptive  Control 

The  uncertainty  in  the  number  of  modes  to  be  selected  for  adequate 
control  of  LSS  raises  some  special  difficulties.  Conventional  adaptive  con¬ 
trol  schemes  involve  controller  structures  of  fixed  order.  In  the  LSS 
context  it  seems  necessary  to  adapt  the  controller  order  as  well  as  its 
parameters.  The  lattice  structure  is  especially  well  suited  for  variable 
order  modeling  and  control  as  was  discussed  in  Friedlander  (1982), 
Sundararajan  and  Montgomery  (1982). 

The  theoretical  development  of  lattice  filters  for  control  purposes  is 
only  at  a  preliminary  stage  of  development.  Its  applications  to  practical 
control  situations  have  been  very  limited.  A  number  of  issues  need  to  be 
resolved  before  the  potential  of  lattice  filters  can  be  fully  realized, 
including: 

(i)  The  development  of  order -recursive  lattice  models  for  plants 
with  both  poles  and  zeros.  The  work  of  Sundararajan  and 
Montgomery  (1982)  was  limited  to  the  all  pole  case.  It  is  ex¬ 
pected  that  pole-zero  models  will  provide  better  fit  to  the  LSS 
problem.  Some  work  has  been  done  in  this  area  (Friedlander, 
1982)  but  questions  remain  regarding  the  simultaneous  determina¬ 
tion  of  the  orders  of  the  poles  and  the  zeros. 


(ii)  The  development  of  direct  lattice  controllers.  The  main  appli¬ 
cation  of  lattice  filters  so  far  has  been  in  the  area  of  system 
identification.  The  identified  model  can  then  be  used  to  design 
a  conventional  controller.  It  is  possible,  however,  to  use  a 
controller  implemented  in  lattice  structure  and  adjust  its  para¬ 
meters  directly,  based  on  the  observed  data.  This  leads  to  an 
adaptive  controller  structure  in  which  both  the  gains  and  the 
order  are  being  adjusted. 

(iii)  Development  of  a  prediction  error  lattice  filter.  The  current 
versions  of  the  adaptive  lattice  are  of  the  residual  type  (i.e., 
the  predictor  coefficients  depend  on  current  as  well  as  past 
data.  In  the  prediction  error  version  they  depend  on  past  data 
only).  The  control  problem  is  naturally  related  to  the  predic¬ 
tion  error  form.  This  development  is  a  straightforward  exten¬ 
sion  of  available  results. 

Finally,  it  is  necessary  to  test  some  candidate  lattice  algorithms  on 
real  and  simulated  data  to  gain  better  insight  into  their  properties  when 
applied  to  very  high  order  plants. 
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APPENDIX  A 

DERIVATION  OF  GENERIC  ADAPTIVE  SYSTEM 


In  this  section  we  show  how  the  nonlinear  system  (see  Figure  2.1) 


e  =  e^  -  c'  X 
z  »  -  Dx 


X  »  Ax  +  Bz' 6 
0  =  Yze 


(A.1) 


is  representative  of  a  large  class  of  adaptive  control  and  identification 
systems.  In  (A.1),  e(t)  and  e^(t)  are  scalars,  x(t)  is  an  n-vector,  0(t) 
and  z(t)  are  p- vectors.  A,  b,  c,  and  D  are  constant  matrices  or  vectors  of 
appropriate  dimensions,  and  Y  is  a  positive  constant.  A  more  thorough  dis¬ 
cussion  on  the  derivation  of  (A.1)  can  be  found  in  Kosut  and  Friedlander 
(1982,  1985),  or  Kosut  and  Johnson  (198M). 


Adaptive  Control 

Consider  the  model  reference  adaptive  control  (MRAC)  system  of 
Narendra,  Lin,  and  Valavani  (1978): 

•  I 

Plant:  X  -Ax  +bu,x(o)eR 

p  p  p  p  p  p 


y-d+c'x  (A. 2a) 

P  P  P 


Reference  model: 


Control : 


X  -  Ax  +  br,  X  (o)  -  0  E  R 
m  mm  m  m 


y  -  c'  X 
m  mm 


u  -  -0'z 


(A. 2b) 

(A. 2c) 


Filter  (Regressor):  z  »  Aj.z  +  ^fy^p  *  ^fu^p  ”  ^fr*"’  =  0  c  (A. 2d) 


Parameter  Update:  6  »  Yze,  e(o)  e  R^,  Y  >  0 

e  “  y  -  y 
■'p  ■'m 


(A.2e) 


The  external  signals  are  the  reference  command  r(t)  and  the  disturbance 
d(t).  Let  be  a  constant  vector  in  a  subset  of  R*^  such  that  when  the 
fixed-gain  (non-adapti ve)  control 


u  =  ~6^z 


(A. 3) 


is  applied  to  the  plant  the  resulting  system  is  stable  and  in  addition  ex¬ 
hibits  acceptable  performance  characteristics.  Any  system  corresponding  to 
any  0^  e  is  referred  to  as  a  tuned  system  and  is  described  by: 


=  A  x^„  +  b  u„,  ,  y„,  =  d  c’x 


■p*  p"p*  'p  'p*  ’  "^p* 


P  P’ 


z.  *  A-z„  +  b_  y  _  +  b„  u  .b_  r 
»  f  *  fy  P*  fu  p*  fr 


=  ’Qiz* 


(A.M) 


e 


-  y 


m 


Define  the  parameter  error  vector 


Combining  (A. 2),  (A.iJ),  and  (A. 5)  and  comparing  the  resulting  system  to 
(A.1)  gives  the  error  state  as 


X 


(A. 6) 


and  system  matrices: 


A  = 


/  A 

-b  01  \ 

I 

'  ’’fy'p 

P  *  ' 

A.  -  b.  0' 
f  f  u  i 

,  b  -  I 

(A, 7a) 


(c^  0),  D  =  (0  Ip) 


(A. 7b) 


Since  the  tuned  system  is  stable  by  definition,  it  follows  that  A  is 
stable,  i.e,,  all  the  eigenvalues  of  A  have  negative  real  parts. 

Identification 


The  algorithms  studied  in  LJung  and  Sodestrom  (1983)  and  Goodwin  and 
Sin  (1984)  will  be  used  as  representatives.  Consider  the  possibly  unstable 
plant 


c’  X  + 
P  P 


(A. 8) 


A  X 
P  P 


b  u 
P  P 


Xp(o) 


e  R 


with  (A  ,  b  ,  c  )  controllable  and  observable,  to  be  identified  by  using  the 
P  P  P 

linear  parametric  model 


X  ■=  A-X  -  y  +  b,  U  ,  X  (o)  E 

m  fm  fy  fup  m 


where  y  is  set  equal  to  either  y  or  y  ,  A.  is  stable,  and  (A„,[b»„  b 

pmf  iiy 

is  controllable.  In  the  case  when  SL  «  p  and  d(t)  =  0,  a  unique  6^  e  B 

ists  such  that  y  -  y  ,  where  y  -  satisfies: 
p  m*  m* 


y  -  =  Olx 

■'  m»  *  rr 


X  „  =  A.x  »  -  b-  y  +  b„  u 
m*  f  m*  fy  fu  p 


y  =  yp  = 


In  practice  5,  >  p  due  to  unmodeled  dynamics  and  d(t)  ^  0  due  to  sensor 
noise  and  disturbances.  In  this  case  6^  e  can  be  chosen  to  minimize 
either  the  equation  error  norm 


J,(e)  -  l|yj.)  -  y.<.)|| 


y  -  y. 


or  the  output  error  norm, 


Jj(e)  -  ||yp(.)  -  y„(.)||j,  . 


where  the  norm  ll-l]  is  defined  as 


I |f(.)| I  -  (  lim  1  [  |f(t)|^  dt)^^^ 


When  =  p  and  ||d(.)||  =  0,  the  equation  error  solution  is  Identical  to  the 
output  error  solution  and  is  given  by  the  unique  0^^  satisfying  (A.  10). 

Hence,  in  this  case,  >=  »  0.  When  H  >  p  or  ||d(.)||  >  0,  the 

solution?  are  different.  In  general,  output  error  provides  some  filtering 
of  the  disturbance. 

The  on-line  or  adaptive  procedures  for  obtaining  output  error  or  equa¬ 
tion  error  estimates  have  the  following  forms.  For  equation  error 
identification; 


Error: 

e  =  y  -  y 
•’  p  ■'m 

•  %  = 

9'x 

m 

Model ; 

=  A  - 

m  f  m 

b-  y  + 

fy^p 

b-  u 
fu  p 

(A. 14) 

Update: 

0  =  Yx  e 
m 

For  output  error  identification: 


Error: 

e  =  y  -  y 
•'  p  ■'m 

•  =  ®’^m 

Model : 

X  =  A_x  - 
m  f  m 

^fy^m  *  '^fu^p 

i 

Update: 

0  =  Yx  e 
m 

Let  6^  denote  a  minimizing  solution  of  either  (A.11)  or  (A. 12)  and  let 
(A. 10)  describe  the  resulting  tuned  model.  Define  the  parameter  error 

A 

e  =  e  -  6,^  (A. 16) 


Combining  (A. 10),  (A. 14),  and  (A. 16)  gives  the  equivalent  equation  error 


e  *  e,  -  0'z  ,  =  y  -  e;z 


(A, 17) 


0  =  Yze 


which  is  a  degenerate  form  of  (A.1),  i.e.,  the  state  is  not  driven  by  the 
non-linear  term  z' 0  as  it  is  in  the  adaptive  case  (A. 7). 

In  the  output  error  case,  combining  (A. 10),  (A. 15),  and  (A. 16)  gives 
the  equivalent  system 


V-.-; 


e  =  e„  -  e'z,  e„  =  y  -  B’x  - 
*  ’  »  •'p  *  m* 


z  =  A.z  -  b.  u  +  b-  z'0 
f  fu  p  fy 


(A. 18) 


0  =  Yze 


which  is  also  a  degenerate  form  of  (A.1).  Note,  however,  that  in  the  output 
error  case,  the  state  equation  is  driven  by  z' 0  whereas  in  the  equation  er¬ 
ror  case  it  is  not.  Hence,  the  output  error  systems  is  much  more  like  the 
adaptive  system  (A. 7)  and  can  exhibit  the  same  type  of  instabilities  and 
parameter  drift  that  have  been  reported  in  several  studies,  e.g.,  Rohrs  et 
al.  (  1982,  1933),  Reidle  and  Kokotovic  (I98i1). 


*'• 


•V  S".  v'- 


APPENDIX  B 


METHODS  OF  AVERAGING  FOR  ADAPTIVE  SYSTEMS* 


Robci'l  L.  Kosut 
Ini  fgi'dled  Syst  ems.  Inc. 
101  University  Avenue 
Palo  Alto,  CA  qii30l 


Absl  racL 


A  sumiriary  of  methods  of  averaging  analysis 
is  ed  for  continuous- time  adaptive 

Systems.  The  averaging  results  of  Riedle  and 
K  'I9.8e)  and  of  Ljung  M9''7)  are  examined 

and  a'.,  sncwn  to  be  closely  related.  Both 
appr,:.,  n.  s  result  in  a  sharp  slabi  l  ily- instability 
bounda''y  which  can  be  tested  in  the  frepuency 
dom-ir.  -r;;!  i  rite-pr  e' ed  as  a  signal  dependent 
positivity  condition. 


' .  1 ntroduc* icn 

The  developed  in  Kosut  and  Anderson 

‘qsS)  shows  that  the  stability  of  adaptive 
sys'tms  in  the  neighborhood  of  ln<  equilibrium 
'.ra j'-ctori es  is  dependent  on  the  stability  of  a 
sysle.m  of  linear  lime-varying  equations.  These 
•■'SS-h'.  ially  are  a  linearization  of  the  adaptive 
System  and  are  r-ftrred  to  here  as  th<  linearized 
adaptive  system. 

In  a  recent  paper  by  Riedle  and  Kokclovic 
;  classical  method  of  averaging  as 
des.-ripr.;;  (jy  Hslc  (19^0)  was  applied  to  the 
linearized  adaptive  system.  The  result  is  a  sharp 
stabi 1 ity- Ihstabl 1 ity  boundary  determined  by  a 
s.gnal  dependent  positivity  condition.  This 
re.iuit  is  si  gni  fi  eanl.  ly  weaker  than  the  SPR 
strictly  pcsilive  real)  ccndilion  required  in  the 
P'oof  of  stability  of  adaptive  systems,  e.g.. 


Narendra,  Lin,  and  Val.ivannl  (  1980),  Landau 
(19V9).  In  order  to  apply  the  averaging  theory  to 
obtain  tnis  result,  the  linearized  system  has 
first  to  be  decoupled  into  slow  (parameter)  stales 
and  fast  stales.  It  is  this  transformation  which 
is  essential  to  the  averaging  approach  and  is  a 
major  contribution  in  the  Riedle-Kokotovi c  method. 

Averaging  has  also  been  applied  to  the 
counter-example  of  Rohrs  et  al .  (1983)  by  Astror. 
(1983,  1989).  In  this  analysis,  by  "freezing"  the 
parameters,  the  parameter  and  state  equations  are 
decoupled  thereby  obtaining  the  asymptotic 
trajectories.  Both  of  these  averaging  analyses 
assume  that  the  system  is  periodic  or  almost 
periodic,  an  assumption  that  can  be  dispensed  with 
by  introducing  the  notion  of  a  local  (moving) 
average,  Kosut,  Anderson,  and  Mareels  (1965).  In 
the  same  reference,  the  averaging  approach  is  also 
shown  to  be  applicable  to  discrete-lime  systems. 

In  Riedle  end  KokolovSc  (1985),  the 
averaging  approach  is  extended  to  nonlinear 
systems  -  and  generalizes  Astrom's  analysis  by 
Introducing  the  integral  mani f cl  d  which  completely 
separates  the  parameter  and  state  equations.  This 
latter  approach  is  valid  for  the  nonlinear 
adaptive  system,  and  not  just  the  linearized  part. 
Similar  results  can  also  be  found  in  Fu,  Bodson, 
and  Sastry  (1985) . 

An  averaging  method  was  also  developed  by 
Ljung  (1977)  for  use  in  discrete-time  recursive 
parameter  estimation.  The  analysis  shows  that  the 
convergence  properties  of  the  estimates  can  be 
determined  from  the  stability  properties  of  a 
related  set  of  ordinary  differential  equations; 
the  method  usually  referred  to  as  the  ODE 
analysl s . 

In  this  paper  we  summarize  the  results 
obtained  by  Riedle  and  Kokotovlc  (1989,  1935)  and 
show  (heur i 3t i ca 1 1 y )  how  they  are  related  to  the 
local  st.'ibillty  analysis  of  Ko.sul  and  Anderson 
('983,  1989)  and  Ihf  ODE  averaging  approarh  of 
Ljung  (1977). 


•  Rese.ir?!'.  Suppor  ted  by  AFOSR  under  Contrael 
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A  !  ,£.1  ;  KfTi.;?-  r.y.t  om 

A  ]  f  ^l^:  it  ]  .s  utilik^ly  i  N  i»  a  truly  g»'n»rir 

i'l.f.V  .vt  error  syslefn  c.m  hr  fcrin»cl  Lo  CMplure  all 
’  fr  M.n  a  cf  ad-jpliv»^  syiift-rni:,  the  adaptive 

^rn  ah^ywn  in  Figure  1  is  offered  as  a  good 
r  •  ; s.  at  i  on  for  in-.*  pui’posen  cf  analysis.  Tht‘ 
syi*-'-”  '-qUel  ions  .-.n-: 


i r.i,  tl  Stat'i  1  i  t  y  and  Pass!  v  >  t  y 


It  is  cf  interest  to  determine  under  what 
conditions  the  adaptive  er’ror  system 
produr-es  bounded  outputs  ^O.e.v.z)  fcr  al  1  bound*^  d 
initial  paramr-ter  t-rrors  u(c)  c  R^.  This  is  whit 
is  m-.-ant  here  by  "global"  stability.  As  it  turns 
Out,  such  a  result  is  possible  to  prove  provid'd 
that: 


m 


I)  c  SPR  with  gradient  (3.1) 

?)  H  (s)  -1/;?  e  SPR  with  least  squares 

(3.2) 
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r  a-f.  -  IrpnirnL  of  (2.1)  can  be  found  in  Kosuf  and 
.-dlar.dt-r  M9Rj,  1985)  and  in  Kosut.  and  Johnson 
i"-.  In(2.1),e(l)eRisa  measured  error 
rna’.  which  drives  the  paraneler  update  (T.ld), 

'  ■.  S'"  is  th.-  regi’essor,  and  O'l)  c  is  the 
"i.T.ct.r  error  between  the  current  estimate  at  t 
1  a  un-"1  parameter  setting  0,  c  R^.  The 
-•■tlitn  of  S,  is  based  on  complete  knowledge  of 
-  'r'.uai  plant  and  di  si  rubances  .  The  system 
"espor.ding  to  this  setting  is  referred  to  as  the 
■e-:  sysi  em.  The  signals  e,(t)  c  R  and  z,(t)  c 
are  cutputs  of  the  tuned  system,  and  are 
r.rred  to  as  the  tuned  error  and  tuned 
g"-sso- ,  respectively.  The  signal  v(t)  c  R  can 
e-;gar  ded  as  the  adaptive  control  error. 


The  iierators  H  and  H  are  dependent  on 
1,  ai.d  describe  how  V  effects  the  error  and 
r- signals.  We  assume  here  that  and 
are  1 ;  near- 1  i  m.  -  i  nvar  i  ant  (LTI)  with  sifable 
proper  transfer  functions  is)  and  (s).  This 
w-,..d  arise,  fcr  example,  when  the  plahJ!  to  be 
:;r.'.ro.led  is  LTI  and  Lh€f  adaptive  controller  is 
iin-je  in  tn,.  adaptive  parameters.  The  stability 
'f  H  ind  H  IS  ■(  consequ.’nce  of  the  definition 
cf  w,’^as  the“{uri.’d  parameter  setting. 

Tr^e  operalor  f  depends  on  the  choice  of 
pi'treler  update  algorithm.  We  will  restrict 
attention  here  to  the  following  representatives: 


(rzlii)  •  cz(t) 


R^-Curs i  ve  Least  Squares 


(rzHt)  .  p(i)zft) 


-  p'’it)  .  z(l)  z(t)' 


3)  z,,z,E  and  either  (3.3) 

a)  e,,^,  t-  LjA  L^  (3.1) 

b )  e, E  and  (3.5 ) 
z  E  PE  (persistently  exciting) 


Parameter  convergence  to  a  constant  in  R'^  or  to  a 
well  defined  subset  in  R^,  requires  that  (3.ii)  be 
strengthened  to: 

E  Lj  A  L^,  z,  E  PE  (3.6) 

The  above  results  can  be  found  in  Kosut  and 
Friedlandcr  (1982,  1985)  and  in  Boyd  and  .Sasiry 
(19811).  Although  of  theoretical  significance, 
they  are  not  feasible  to  obtain  in  practice.  In 
the  first  place,  due  to  unmcdeled  dynamics  (Rohrs 
el  al .  1982,  1989),  (s)  e  SPR  is  praticaliy 

impossible  to  achieve^Tn  adaptive  feedback  and 
even  in  some  output  error  identifica-  lion.  (This 
is  not  the  case  in  equation  error  identification.) 
Secondly,  when  e,  e  as  in  (3.5),  it  is  required 
that  z  c  PE  which  can  not  be  guaranteed  in  advance 
since  z  is  i nsi de  the  adaptive  locp.  Case  (3.6) 
which  rE--quries  z,  c  PE  which  is  feasible  to 
establish  -  conflicts  with  e,  ,6.  e  L^a  The 

latter  implies  e,(t)  -»  0  which  can  only  occur  for 
z*  c  PE  -  and  where  there  are  no  unmodeled 
dynamics  which  we  argue  is  net  possible. 

With  these  impossible  to  satisfy  theoretical 
requirements,  it  is  doubtful  that  a  global 
stability  theory  can  be  attained  which  relies  on 
passivity,  l.e.,  condition  (3.1,?).  On  the 
practical  side,  however,  there  is  substantial 
evidence  of  well  engineered  algorithms  that  work 
without  SPR,  e.g.,  Astrom  (1983).  These  do  not 
work  for  all  0(o)  and  for  all  e, ,z,  in  L^,  but 
rather,  for  restricted  magnitudes  and  signal 
spectruras.  For  example.  If  H  (s)  is  SPR  for  ui  £ 
UL  then  it  is  expected  that  file  adaptive  system 
will  be  well  behaved  provided  there  is 
insignificant  excitation  above 


P(0:  -  Pi0''>0 
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Py  restrict,ing  the  rndgnit-udf  of  0(o)  and  the 
rajgnit  ude  and  apt-clrum  of  z,  (1)  and  e,(t.),  it,  la 
possible  to  obtain  conditions  to  prove  local 
stability,  e.g.,  Kosut  and  Anderson  (19fl3,  1985). 
The  local  stability  property  hinges  on  two 
premises:  M)  the  error  system  trajectories  are  in 
a  (not  necessary  small)  neighborhood  of  the  tuned 
solution,  and  (?)  the  linear  time  varying  system 
winch  maps  w  -•  0  as  given  by 


9  -  '  (rz,)(t)  (z,{.)0{.))  e  (rw)(t)  (H.1) 

w  gy  w 

is  L^-  stable,  i.e,,  there  exists  constants  k  and 
b  s.l.  k||w||_^f  b.  The  choi  ce  of  f  comes 

from  ((’.D  or  (2.3)  and  z,  (t)  is  the  tuned 
regressor.  We  can  regard  (A.l)  as  a  linearization 
of  the  update  algorithm.  There  are  several  ways 
to  establish  the  stability  of  (A.l). 

G-adicnt  Algorithm 

We  first  consider  the  case  when  T  represents 
the  gradient  algorithm,  i.e.,  (rz,)(t)  «  cZj,(t) 
with  t  >  0. 

In  Anderson  (1977),  it  is  shown  that  if  H 
(s)  e  SPR  and  z,  c  PE,  then  for  all  e  >  o,  w  -•  6'^ 
is  exponentially  stable,  and  hence,  L^-  stable. 

In  Anderson  ct  al .  (1935),  if  H  (s)  •  H  (s)  ♦ 

ev  ev 

‘sAiSl,  H^,^(s)  c  SPR,  a{s)  is  stable,  and  z,  e  PE 
then  for  sufficiently  small  e  and  |ji,  |l,.  w  •*  0 
is  still  exponentially  stable,  and  hence  is  bi¬ 
stable.  This  latter  method  relies  on  loop- 
transformations  and  applicalton  of  small  gain 
theory. 


B.  Recursive  Least  .Squares  Algoriltim 


In  this  case  we  have  from  (2.3)  that 
(rz,)(t)  -  P,(t)z,(t) 

^  P,(t)’'  -  z,(t)z,(t)',  P,(0)>0. 
When  z,  £  PE  there  exists  a  >  0  such  that 


P,(t)  '  -  P,(o)  '  *  /  z, ( T)t^( T ) ' di 

0 

>  at  .  1 


Thus,  it  is  convenient  to  define 
R(t)  -  P.(t)"’  for  t  >  o 

t  * 


R(t)”’  .  t  P.(t)  <  -  1 
*  £1 


and  we  can  write  (A.l)  and  (2.3)  as, 


0  -  1  R'Vw  z.H^^(z.O)] 


R  -  -  (2),Z*  "  R) 


When  H  (s)  -  is  not  SPR  we  can  now  follow  Ljung 
(1977)®and  for^t  >  s  and  s  sufficiently  large, 
replace  the  rlght~hand  side  by  its  average. 

Letting  "overbar"  denote  average  (assuming  it 
exists)  we  have: 


b  % 

*  ^  *  •  •  •  • 
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G(t)  .  R  (w  -  (2,H^^2;)0) 
R(t)  -  -J-  (z,zi  -  R) 


Another  approach  is  to  use  averaging.  In 
Riedle  and  Kokotovic  (198ii)  it  is  shown  that  if  z, 
e  PE  with  th“  Fourier  series  representation 


z,  f  t )  ■  I  Q  ( )  e 


and  if  the  eigenvalues  of  the  real  matrix 

B  -  I  a(u.^)  o(-W|^)'  f^ev^'-’^k^ 
k 

all  have  positive  real  parts,  then  for  ail 
sufficiently  small  t  >  o,  w  -►  0  is  exponentially 
stable,  and  hence,  Li-stable.  Moreover,  if  any 
one  eigenvalue  of  B  has  a  negative  real  part,  then 
w  -  0  is  exponentially  unstable.  Hence,  there 
exists  w  £  Li  s.t.  |Q(L)|-*  s’  as  t  " 
exponentially  fast.  It  is  obvious  then  when 
H  (s)  is  not  SPR,  but  only  approximately  so,  then 
t^e  Riedle  Kokotovic  result  provides  a  sharp 
5tabil i ty- 1 nstabi 1 i ty  boundary.  Note  than  when 
(si  is  SPR  and  z,  c  PE  we  have  from  Anderson 
(T§77)  that  w  ->  0  is  exp.  stable  for  all  t  >  o. 

On  the  other  hand,  the  result  in  Anderson  et.  al. 
'1935)  remains  valid  for  H  (s)  t  SPR  (a(s)  •  o) 
because  then  c  >  o  is  bounSed  above  by  infinity. 


Integrating  from  s  to  s  •*  T,  T  >  o,  gives 


.  -1  _ 

[0(3*T)-O(s)]/;  dt/t  •  R  (w-  (z.H  zllO) 

*  ev 

5 


[r(3*t)-r(s)]/;  dt/t  -  -  R  (11.7) 


Now  change  time  scales  s+T-»T+ii,  it-/ 

dt/t  and  letting  s  -*  •  gives  the  differential 
equations: 


Oj^(t)  -  R^(t)'''[w  -  B  e^(T)] 


R  (t)  -  z.z;  -  R  (i) 
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with  ttic-  gradient  algorivtim  (?,?),  lei 


with  B  -  given  by  (‘i.3).  These 

eqmllons  des'^ribe  ihe  asymptclic  behavior  of 
(tt.l)  in  Just  the  same  way  as  they  do  for 
di  3crele"lim»'  (Ljung  19^7). 

Observe  that  in  (9.9)  as  i  ■*  ".  >  )  ■* 

e»^'i  •  Thus,  when  H^^(s)  ~  -  is  hol  SPil  and  t 

PE  wi Ih  Ecurier  representallon  (9.P)  the 
asymptotes  are  st.able  if 


♦  d)  -  0(1) 
Y(t)  -  c 


(B.5) 


and 


f(t,0,x)  -  2,{t)e,(L)-  0,(1)  X  ♦  c'x  D  X 


Re  1  (L)  >  o 


(M. to) 


Q,(t)  -  z,(t)c'  ♦  e,(t)D 


(5.6) 


g(t,0,x)  -  b(2,(l)  -  Dx)'0 

where 


L  -  (z,zj)  '  B 

-  (I  a(  u.^)a(  "ui^) ' )  '  I  a(  u)^)a(  Uj^) 'H^^(-Ju^) 

l<  ^  (if.tt) 


With  the  recursive  least  squares  algorltlim  (2.3). 
define: 

R(l)  -  I  Pd)'’  (5.7) 


If  Re  H(jie  )  a  p  >  0  at  low  frequencies,  and  if 
|a(u  )|  is'^small  at  frequencies  whore  Re  H(JU|^)  < 
0,  tJien  Re  1(L)  *  p.  Thus,  all  parameter 
asymptotes  have  a  uniform  rate  of  convergence 
which  is  not  the  case  for  the  gradient  algorithm 
with  a  time-invariant  gain. 


and  let 


where  the  oprerator  collR)  stacks  up  the  columns 
of  the  matrix  R  to  form  a  vector.  Thus, 


5 •  Av-raginp:  A  More  General  Approach 


In  this  section  we  will  establish  a  general 
form  of  the  adaptive  error  system  (2.1,3)  which  is 
useful  for  application  of  averaging  methods.  The 
first  step  is  to  transform  (2.1,3)  Into  a  set  of 
nonlinear  time-varying  differential  equations.  To 
do  this  observe  that  if  H  (s)  are  strictly  proper 
functions  (a  convenient  illfustrative,  but  not 
necessary,  assumption)  then  we  can  write 


H  (s )  •  o' ( si -A )  ' b 
ev 

H  (s)  -  D(sl-A)''b 


(5.1) 


where  A  £  r"’'"',  b  c  r",  c  c  r",  D  c  rP’'",  with 
;A,b.[c  D'])  3  minimal  representation.  Also,  Re 
i(A)  <  0  reflecting  the  fact  that  and 
are  stable.  The  error  system  (2.1)  is  then 
equivalently  expressed  as 


e 


e. 


c'x 


z  -  z,  '  Dx 
X  -  Ax  ■»  bz'3 


(5.2) 


0  -  (r7,)e 


Sy  eliminating  the  variables  e  and  z  we  can  reduce 
:5.2:  to  the  coupled  state-space  description: 


(5.3) 


f(t,4i,x)  • 

R  ’(z,(l)e,(t)  -  Q,(t)  X  ♦  c'xDx) 

col(z,(t)z,(t)’  -  z,(t)(Dx)’  -  Dxz,(t)'l 
♦  Dx(Dx)'  -  R) 


(5.9) 


g(t,*,x)  •  b(z,d)  -  Dx)'0 


(5.10) 


The  coll.)  operator  was  used  by  tjung  (197S)  to 
develop  the  discrete-time  version  of  (5.3).  (5.9). 


Heuristics:  The  Integral  Manifold 

The  basic  idea  in  the  application  of 
averaging  methods  to  (5.3.9)  is  to  see  what 
happens  when  7(t)  is  small.  Essenti allyl .  ♦(!) 
slows  down  and  we  can  repalee  the  right  hand  side 
of  (5.3)  with  its  average,  l.e.. 


^  -  Y(t)f((f)  (5.11) 

1  T 

where  fl^)  •  Hm  ^  ^  r(t  ,x(  t  ,♦ )  )dt  (5.12) 

T-*“  0 

assuming  the  limit  exists,  (Such  is  the  case,  for 
example,  when  f(t,4i,x)  and  gd.if.x)  are  periodic 
in  t  for  all  bounded  ♦  and  x).  To  arrive  at 
(5.11)  formally  requires  the  Introduction  of  the 
integral  manifold  as  suggested  by  Reldle  and 
Kokotovlc  (1985)  [see  Hale  (1969)  for  discussion 
of  the  Integral  manifold] 

The  integral  manifold  M  of  (5.3,9)  is  the 

set. 


i  .  7{t)  f(t,((l,x) 
X  -  Ax  ■*  g(t,(ti,x) 


(5.9) 


M  -  (t,»,x  :  x(t)  -  h(l,(Kt)) 


(5. '3) 


Theorem 


dy  subs!  i  I  ul  1  rig  x  »  h(l,if)  Into  l-be 

imnifold  f  unfit  on  hd.iti)  Is  seen  lo  satisfy  the 
partial  di  f  f  orc-nt  i  al  equation 

-  >;t)  rit.D.h)  -  Ah  <  g{t.(ti,h)  (S.IA) 


Let  denote  3  solution  of 

r(»°)  -  o 


Whenever  y(t)  is  sufficiently  small,  a  reasonable 
approximation  is  to  hCt.iti)  Is  given  by  h^(t,(|i) 
wiilch  is  the  solution  to 

_ o  -  Ah^  *  gCt.if.h^) 


-  F(u)h  ♦  bz,(t)'G 

o  • 


(5.15) 


where  the  last  line  follows  from  (5.6)  with 


FiO)  -  A  -  bO'D 


(5.16) 


In  (S.ltt),  Q  and  t  are  regarded  as 
independent  variables  and,  hence,  we  can  define 
the  stabilizing  parameter  set 

D  -  10  c  H^r  Re  1(F(0)  <  o)  (5.17) 

3 

Thus,  for  Y'lj  sufficiently  small,  we  can  refer  to 
hit,!})  with  G  c  D  as  the  stable  manifold,  which 
we  will  approximate  by  h^(t,G),  G  e  D^. 

The  final  transformation  on  (S-B.*!)  is 
obtained  by  examining  the  behavior  of  (iji.x)  In  the 
neighborhood  of  the  stable  manifold.  Introduce 
the  error  state. 


C  -  X  -  h(t,») 


(5.18) 


Using  (5.18),  and  (S.B,**).  with  the  approximation 
h^(t,G)  for  h(t,<)i),  we  have 

<p  ■  >(t)f(t,p,hg(t,Q)e5;)  (5.19) 


and  define  the  matrix. 


Provided  Re  1(C)  W  o; 

(1)  If  Re  KG)  <  o  then  Is  u.a.s. 

(il)  If  max  Re  1(C)  >  o  then  Is  unstable. 


Application  to  Gradient  Algorithm 

Applying  this  result  to  (5.5,6)  with  the 

gradient  algorithm  and  with  z„  c  PE  and  z,e, 

-  o,  gives  C  «  -B  with  B  from  (I4.3).  Since  ir(t)  - 
E  >  o.  we  can  only  conclude  that  If  Re  KB)  >0 
and  £  Is  sufficiently  small,  then  0(t)  approaches 
an  £-neighborhood  of  G  -  o  as  t  -►  ».  Provided  of 
course  that  0(t)  c  long  enough  for  transients 
to  die  out,  which  is  unprovable  as  yet  in  general. 


Application  to  Recursive  Least  Squares  Algorithm 

Under  ths  same  conditions  and  with  the  same 
provisions  as  above,  C  --L  with  L  from 

(9.11).  This  time,  since  Y(t)  •  1/t  o  as  t  ■»  “, 
we  can  conclude  that  If  Re  1  (L)  >  o,  then  G(t)  ■» 
o  as  I  ■»  «  at  a  rate  1  /t. 
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APPENDIX  C 


Fundamentals  of  Control  for  Large  Space  Structures: 
Robust  and  Adaptive  Design 


Abstract 

This  report  gives  a  summary  overview  of 
robust  and  adaptive  control  design  for  LSS. 

It  is  one  section  of  a  larger  internal  report  in 
progress. 

11-2.2.7.4.1  Control  Design  for  HAC/LAC 
Architecture.  In  this  section  we  will  discuss  the 
steps  involved  in  control  design  for  the  HAC/LAC 
architecture.  Although  the  architecture  is  special- 
ired.  the  control  design  methodology  is  not  and  can 
be  quite  general.  We  will  discuss  three  methodolo¬ 
gies  for  design:  (1)  an  LQG  based  methodology 
whose  genesis  is  the  ACOSS/VCOSS  programs, 
and  (2)  a  more  recent  approach  involving  what 
is  known  as  "Q-parameterization  and  “’Hoe- 
optimization".  These  latter  methods  are  frequency- 
domain  oriented  rather  than  state-space  oriented 
like  the  LQG  approach.  {31  We  will  also  discuss  an 
adaptive  control  strategy  which  can  be  utilized  for 
online  self-tuning.  We  refer  to  this  approach  as 
"adaptive  calibration".  This  approach  has  been 
developed  by  ISI  under  an  on-going  research  con¬ 
tract  with  the  Directorate  of  Aerospace  Sciences 
in  AFOSR 

Limitations  of  Design.  Independent  of  the 
design  method,  the  defining  characteristic  of  the 
vibration  control  problem  is  there  are  an  infinite 
number  (theoretically)  of  elastic  modes,  with  low 
natural  damping,  and  the  controller  bandwidth 
extends  over  a  significant-number  of  these  modes 
{Fig  II. 2-28)  The  low  frequency  modes  interact  not 
only  with  the  attitude  controller  but  contribute 
directly  to  the  deformation  geometry  of  the  structure 
which  itself  may  require  accurate  control  Proper 
control  synthesis  requires  that  performance  criteria 
be  precisely  formulated  or  the  control  problem  is 
ill- posed. 


Fig.  11.2-28  Flexible  Structure  Mode  Location 
and  Controller  Bandwidth 

The  control  design  approach  must  properly  han¬ 
dle  the  poorly  known  higher-frequency  modes  by 
not  destabilizing  them  while  controlling  the  low- 
frequency  modes.  Indeed,  no  matter  where  the 
controller  roll-off  frequency  is  situated,  the  infinite 
nature  of  the  modal  spectrum  implies  that  there  will 
be  modes  within  and  beyond  the  roll-off  region,  fur¬ 
thermore.  destabilization  is  likely  and  almost  certain 
to  occur  in  the  roll-off  region,  a  situation  which  can 
only  worsen  for  closely  packed  modes  and  low 
natural  damping  This  phenomenon  sometimes 
referred  to  as  “spillover"  is  one  of  the  most  crucial 
problems  faced  by  the  control  designer.  In  more 
general  terms,  spillover  can  be  viewed  as  an  aspect 
of  the  problem  of  robust  control  design:  this  will  be 
discussed  more  in  a  later  section. 

Modeling  of  Flexible  Spacecraft.  A  central  issue 
in  the  active  control  of  space  structures  is  the  devel¬ 
opment  of  "correct"  mathematical  models  for  the 
open-  and  closed-loop  dynamical  plants  Programs 
such  as  NASTRAN  and  SPAR  are  the  primary  cur¬ 
rent  tools  foi  generating  dynamical  models  of  con¬ 
ceptual  spacecraft  whose  structure  cannot  be  ideal¬ 
ized  by  simple  models  of  beams,  plates,  and  beams 
with  lumped  masses 

Finite-element  structural  programs  generally  pro¬ 
vide  the  control  designers  with  a  set  of  modal  fre¬ 
quencies  and  a  set  of  mode  shapes  (eigenvectors) 
corresponding  to  appropriate  boundary  values 
(e  g  .  free-free  modes)  These  eigenvectors  are 


given  in  discretized  form,  i  e  .  a  set  of  modal  dis¬ 
placements  in  the  X.  y  and  z  directions  at  each 
nodal  station  In  some  cases,  modal  rotations  are 
also  required  In  addition  coordinates  and  a  "map" 
of  the  structure  s  nodes  must  be  provided  to  allow 
the  reconstruction  of  physical  displacements  in 
terms  of  their  modal  expansions 

The  important  point  here  is  that,  for  any  nontri¬ 
vial  flexible  satellite  configuration,  the  volume  of 
information  is  so  large  that  the  data  handling  must 
remain  entirely  within  the  computer  and  its  mass- 
storage  facilities  Development  of  this  database,  in  a 
form  usable  by  control  synthesis  software,  is  a  fun¬ 
damental  necessity  for  the  synthesis  and  evaluation 
of  complex  controls  which  require  modal  trunca¬ 
tion.  actuator/ sensor  location  and  type  changes, 
and  evaluation  of  system  performance  for  parame¬ 
ter  and  system  order  changes.  Preparation  of  a 
structure  for  controls  is  a  major  part  of  the  overall 
effort  required  to  develop  structural  control  systems. 

Nonlinear  Models.  For  single-body  monolithic 
structures,  the  fine-pointing  attitude  dynamics  arc 
subsumed  in  the  rotational  rigid  body  modes 
included  in  the  modal  matrix.  When  only  "small" 
motions  of  a  space  structure  are  being  considered, 
the  conventional  linear  structural  dynamics  analyses 
(NASTRAN  and  SPAR)  are  adequate,  and  the 
rigid-body  modes  are  formally  handled  together 
with  the  elastic  modes,  even  though  the  actuators 
necessary  to  control  them  will  be  different,  in  gen¬ 
eral.  from  those  used  to  control  elastic  vibrations. 
When  larger  attitude  angles  need  to  be  considered, 
if  the  angular  rates  remain  small,  the  linear  equa¬ 
tions  are  still  applicable  provided  that  the  rigid  body 
modes  are  now  given  in  terms  of  three  attitude 
angles  which  then  constitute  the  first  three  modal 
coordinates  The  displacements  are  then  interpreted 

as  the  linear  deformations  of  the  structure  with 
respect  to  the  rotated  frame  This  procedure 
removes  the  kinematic  nonlinearitics  resulting  from 
the  linear  stretching  of  the  structure  under  the  classi 
cal  rigid-body  modes  However,  for  large  angular 
rates,  nonlinear  dynamic  effects  have  to  be  mod 
eled.  even  though  structural  deformations  can  still 
be  represented  by  linear  equations 
Two-Level  Control  Design:  The  HAC/LAC 
Methodology.  The  two- level  approach  consists  of 
a  wide-band.  low-authority  control  (LAC)  and  a 
narrow-band,  high-authority  control  (HAC).  HAC 
provides  high  damping  or  mode-shape  adjustment 
in  a  selected  number  of  modes  to  meet  performance 
requirements.  LAC.  on  the  other  hand,  introduces 


low  damping  in  a  wide  range  of  modes  for  maxi¬ 
mum  robustness.  Figure  II  2  29  shows  the  control- 
design  procedure  with  integrated  LAC  and  HAC 
designs 


Fig.  11.2-29  Analytical  Control-Design 
Procedure 


LA(  !<>  usuallv  miplomcnicd  with  colocated  sen- 
sot<>  and  actuators  However,  the  theory,  based  on 
the  work  of  Auhruti  is  applicable  to  multiple 
actuators  sensors  with  cross  feedback  and  possible 
filters  |Auh  1] 

HAC  uses  a  collection  of  sensors  and  actuators 
not  necessarily  colocated  Selecting  the  increase  in 
damping  ratio  is  realized  by  any  number  of  meth¬ 
ods,  including:  LQG  w'ith  frequency  shaping. 
Q-parametization .  or  Hoo-optimization.  These 
methods  provide  roll-off  over  desired  frequency 
regions  HAC  may  destabilize  modes  not  used  in 
the  design  LAC  is.  therefore,  necessary  to  "clean 
up"  problems  created  by  HAC. 

The  need  to  integrate  HAC  with  LAC  is  shown  in 
Fig  11.2-30.  HAC  is  based  on  models  valid  over  a 
limited  frequency  region.  It  produces  large  increases 
in  damping  ratio  and  disturbance  rejection  in  the  fre¬ 
quency  range  of  interest.  The  effect  of  the  HAC 
controller  on  modes  not  used  in  the  control  design 
and  outside  the  controller  bandwidth  may  be  stabi¬ 
lizing  or  destabilizing.  LAC  is  designed  to  provide 
protection  such  that  adequate  damping  is  provided 
in  the  mode  most  adversely  perturbed  by  HAC. 
With  reference  to  Fig.  II. 2-30.  the  LAC  moves  the 
entire  uncertainty  region  above  the  zero  level  damp¬ 
ing  ratio. 

In  the  next  few  sections,  a  more  in-depth  discus¬ 
sion  of  the  blocks  in  Fig  II. 2-29  will  be  presented,  in 
particular,  actuator/ sensor  location,  model  and 
controller  reduction  methods,  and  HAC/LAC  syn¬ 
thesis  These  methodologies  rely  on  certain  proper¬ 
ties  of  feedback  control:  this  raises  the  issue  of 
robust  control  design  which  is  fundamental  to  the 
whole  design  philosophy  of  feedback,  especially  for 
LSS.  and  this  will  be  discussed  first 
Robust  Control  Design.  This  section  will  describe 
how  to  evaluate  the  robustness  of  a  control  design. 
The  evaluation  is  independent  of  the  methodology 
used  to  achieve  a  particular  design  To  illustrate  the 
technique  we  will  consider  the  robust  control  prob¬ 
lem  of  vibration  suppression  with  unmodcled  high- 
frequency  dynamics  Figure  II  2-31  shows  the  con¬ 
trol  system  where  P(s)  is  the  plant  transfer  function 
matrix  from  actuator  inputs  to  LOS  sensor  measure¬ 
ments  and  where  C(s)  is  the  controller  transfer 


HAC  MODI  I 
eAHDMlOTH 


Fig.  II.2-3D  heed  to  Integrate  High-Authority 
Controller  (HAC)  and  Low- 
Authority  Controller  (LAC) 
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Fig.  11.2-3 1  Vibration  Suppression  Control 
System 

function  matrix.  Neglecting  the  rigid  body  modes  in 
P(s)  and  assuming  infinite  bandwidth  sensors  and 
actuators. 

P(*)  r  1  C.  Is) 
k  =  i 

where 


Suppose  that  n  of  the  modes  are  known.  Let  P^ls) 
denote  the  known  part  of  P(s). 

For  example.  Pplsl  can  be  obtained  from  P(sl  by 
modal  truncation,  i.e..  the  first  n-modes  of  P(sl  are 
retained.  One  can  ask  the  question:  is  this  the  best 
choice  foi  a  given  model  order  n?  In  general,  it 
depends  on  what  is  meant  by  "best."  For  closed- 


loop  contrt>l  It  IS  usually  better  to  retain  those 
n  rnocJes  which  most  affect  the  closecf-loop  perfor¬ 
mance  How  tcj  select  these  modes  will  be  discussed 
in  the  section  on  model  reduction 

Assuming  the  modes  have  been  selected,  define 
model  error  as 

.its)  =  P(s)  P  (s)  l  C.  (s) 
kc:; 

n 

observe  that  6(s)  is  stable  because  both  P(s)  and 
P,,(s)  are  stable.  Hence,  it  can  be  shown  that  the 
closed-loop  system  is  stable  if. 

where  Q,(s)  is  given  by  Qr,(s)  =  C(s) 
[I  P„(s)C(s)]  ■  '  and  where  "of)  denotes  the  maxi¬ 
mum  singular  value  of  the  matrix  argument.  The 
quantity  is  referred  to  as  the  “stability  mar¬ 

gin".  hence,  the  subscripts  “sm".  (See  (Doy.  IJ. 
[Kos  1]) 

The  stability  robustness  test  depends  on  the  loca¬ 
tion  of  uncertainty.  Additive  perturbations  such  as 
those  just  discussed  result  in  the  test  as  shown.  The 
table  in  Fig,  11.2-32  shows  a  variety  of  stability  mar¬ 
gins  corresponding  to  generic  forms  of  model  error. 
In  Fig  11.2-32.  P  =  plant.  C  =  control.  M  nominal 
model,  and  ^  =  model  error.  The  stability  margin  is 
expressed  as  a  function  of  C  and  M  which  are 
known  quantities  Examples  of  some  model  error 
tests  are  shown  in  Fig.  11.2-33  for  the  CSDL  #2 
VCOSS  model 

Performance  Robustness.  The  stability  robust¬ 
ness  tests  can  be  extended  to  evaluate  performance 
robustness  to  model  error  The  evaluation  is  deter¬ 
mined  by  how  performance  is  measured.  Consider 
the  closed- loop  system 

ytt)  =  His)d(t) 

with  H(s)  the  closed-loop  transfer  function. 
Althoi  gh  d(t)  is  not  precisely  known,  it  can  be  con¬ 
sidered  as  the  output  of  a  weighting  filter  W(s) 
driven  by  “noise"  w(t),  then  d(t)  =  W(s)  w(t). 

Output  performance  versus  input  is  shown  in 
Fig  II  2  34  The  table  is  derived  by  boundtng  the 
“noise"  input  wft)  in  terms  of  power,  energy,  and 


magnitude  Except  for  the  mag  in-mag-out  bound, 
all  other  bounds  depend  directly  on  the  frequency 
dependent  quantity  IHljwlWlju.'))  A  natural  fre 
quency  domain  performance  criterion  is  then 

a[HI)  |iW(Mj  *  PU) 

where  p(u.-|  is  selected  on  the  basis  of  power, 
energy,  and  magnitude  specifications  on  the  output 
signals  In  terms  of  model  error  performance  specifi¬ 
cation  IS  satisfied  if 

0  [pftuill  ■  t,  (-) 

L  J  PM 

with  6pn,(u>)  is  the  performance  margin  given  by 

and  where  is  the  performance  of  the  nominal 
closed-loop  system  H^fs)  with  no  model  error. 
Then. 


which  must  always  be  smaller  than  p{u>)  in  order  for 
6p„(w)  to  be  meaningful.  Note  that  6p„(w) 
as  would  be  expected  since  performance  includes 
stability  As  before,  the  location  of  uncertainty  mod¬ 
ifies  the  calculation  of  Spn-lw). 

Usefulness  of  Stability /Performance 
Robustness  Tests 

The  stability  'performance  robustness  tests  are 
indispensible  in  obtaining  a  realistic  preliminary 
design.  They  arc  used  in  a  number  of  places  in  the 
design  cycle  to  establish  the  HAC  ^LAC  gains,  effect 
of  actuator /sensor  dynamics,  and  the  criteria  for 
mode!  and  controller  reduction,  which  will  be  dis¬ 
cussed  in  the  next  section  The  tests  arc  also  invalu¬ 
able  in  establishing  criteria  for  online  system  identift- 
cation  and  control,  which  will  be  discussed  later  on 
in  this  section . 

Model  Reduction.  In  general,  the  requirements 
for  model  reduction  for  active  control  of  large  space 
structures  must  include  the  following: 

1.  The  reduced  model  should  be  suitable  for  con¬ 
trol  design  and  synthesis  It  should  incorporate 
all  features  critical  for  the  selection  of  a  feedback 
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structure  and  control  gains. 

2.  The  reduced  model  should  accurately  incorpo¬ 
rate  actuator  effectiveness,  sensor  measure¬ 
ments  and  disturbance  distribution  [ACOSS.  1). 

3.  The  dynamical  characteristics  of  interest  in  the 
structure  should  be  represented  in  the  reduced 
model. 

A  basic  methodology  for  model  reduction  which  has 
been  used  successfully  in  ACOSS,  VCOSS  and  a 
number  of  other  programs,  internal  balancing,  is 
now  described.  Other  approaches  also  exist  which 
will  be  examined  in  the  program 
Internal  Balancing 

To  determine  the  most  important  modes  for  con¬ 
trol  design,  many  criteria  must  be  considered  includ¬ 
ing  controllability,  disturbability.  observability  in  per¬ 
formance.  and  observability  in  the  measurements. 
Any  mode  which  is  highly  controllable,  observable, 
and  disturbable  must  clearly  be  included  in  the 
design  model,  however  highly  controllable-but- 
unobservable  modes,  for  example,  are  difficult  to 
judge  Moore  (Moo  1]  has  developed  an  "internal 
balancing"  approach  whereby  asymptotically  stable 
linear  models  are  transformed  to  an  essentially 


unique  coordinate  representation  for  which  control¬ 
lability  and  observability  rankings  arc  identical.  The 
definition  of  internally  balanced  coordinates  follows: 

Pet.  An  asymptotically  stable  model 
I  X  =  Ax  ^  BuJ  is  internally 
jy  =  Cx  I 

balanced  over  lO.'^J  IFF 


At 

Ce  dt  =2, 


■  I  1 

where  ^  diag  .(r^c^. ..a^  i 
^  ^12  n 


I  >  i.  >  0^ 

'  '  I 

Notice  that  the  balanced  representation  is  such 
that  the  controllability  Gtamian  and  observabilitL 
Gramian  are  equal  and  diagonal.  The  o,'s  are 
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termed  "second-order  modes,”  In  general,  the 
required  transformation  "scrambles”  the  original 
coordinate  system  such  that  the  physical  meaning  of 
the  states  is  lost. 

However,  for  lightly  damped  structural  models 
with  decoupled  dynamics,  the  internally  balanced 
coordinate  representation  is  approximately  equal  to 
a  scaled  representation  of  the  model  states.  Thus  it 
is  possible  to  write  approximate  formulae  for  the 
states  in  terms  of  the  original  model.  Three  modal 
rankings  are  considered: 

•  disturbance  inputs  to  LOS 

•  actuator  inputs  to  LOS 

•  actuator  inputs  to  sensor  outputs 

These  "second-order  modes"  rankings  give 
important  evaluations  about  which  modes  to  retain 
and  validity  of  actuator/ sensor  placement.  These 
rankings  are  shown  in  Fig  II. 2-35  along  with  LOS 
modal  cost  |Gre.  1)  computed  using  the  colored 
noise  disturbance 

Here  the  absolute  values  of  the  modal  costs  (for 
the  VCOSS  1  model)  are  used  The  RMS  second- 
order  modes  and  modal  costs  are  plotted  versus 
mode  number  in  Fig,  11.2-35.  Immediately  ident 
is  the  clustering  of  these  modal  phenomena.  The 
disturbance  effect  as  seen  through  the  line-of-sight  is 
constrained  to  clusters  of  modes  as  is  the  ability  to 
measure  and  control  the  model.  The  coincidence  of 
the  controllable  clusters  and  disturbable  clusters 
indicates  a  favorable  actuator/'sensor  configuration 
for  the  problem 

Frequency  Weighted  Balanced  Realizations. 

Balanced  realization  model  reduction  can  be 
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Fig.  11.2-35  Open-Loop  Model  Analysis 

extended  to  finding  a  reduced  model  P^fs)  of  a  high 
order  model  P(s)  such  that 


)  <  1 


where  W^fs)  and  W,(s)  are  output  and  input  fre¬ 
quency  dependent  weighting  matrices.  These  can 
be  chosen  to  reflect  closed-loop  requirements  on 
model  error,  vis  a  vis.  frequency-domain  stability 
and  performance  margins  For  example,  stability  of 
the  closed-loop  system  with  C(s)  designed  from 
P„(s)  is  guaranteed  if: 

Wjs)  =  I 

W,(s)  =  W,(s)  =  C„|s)(UP,(s)C„(s))-i 
The  problem  is  that  Wjs)  is  dependent  on  P„(s) 
which  is  unknown  The  let  out  is  that  its  shape  is 
partially  determined  by  the  performance  spec;  thus, 
we  can  make  an  initial  guess  This  technique  is 
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rcfoticd  to  as  '  advanct-d  loop  shaping  "  This 
involves  an  iterative  problem  which  is  solvable  via 
successive  approximation 

Compensator  Order-Reduction.  An  alternate  to 
plant  order  reduction  is  to  design  a  high  order  com¬ 
pensator  and  then  reduce  the  compensator  order 
Let  C(s)  denote  a  high  order  compensator  of  order 
N  designed  to  control  Pis)  of  order  N  or  larger  Let 
C^is)  denote  a  reduced  version  of  C(s)  to  order 
n<N  Motivated  by  the  stability  robustness  theory, 
view  C(s)  -  C„(s)  as  a  perturbation  Hence,  the 
closed-loop  system  with  P(s)  and  C^is)  is  stable  if: 

sup  5  j  vv(i^)  [cti..)  -  j  <  1 


where 

W(S)  =  (1  -r  Pis)  Cis))-  '  Pis) 

The  weight  W(s)  is  stable  because  the  high  order 
control  Cis)  stabilizes  the  closed  loop  system.  In  this 
case  Wis)  is  known  and  we  can  apply  internal  bal¬ 
ancing  to  find  Cpis)  The  disadvantage  to  this 
method  is  that  it  is  necessary  to  find  a  high-order 
compensator  The  advantage  is  that  once  it  is 
found,  internal  balancing  applies  immediately  since 
the  weights  are  known  On  the  other  hand,  direct 
plant  order  reduction  does  not  involve  control 
design  for  the  high  order  plant,  but  does  involve  an 
iterative  process  since  the  weights  are  functions  of 
the  (unknown)  reduced  model. 

Low-Authority  Control  Design.  Low-authority 
control  (LAC)  systems,  when  applied  to  structures, 
arc  vibration  control  systems  consisting  of  distrib¬ 
uted  sensors  and  actuators  with  limited  damping 
authority  The  control  system  is  allowed  to  modify 
only  moderately  the  natural  modes  and  frequencies 
of  the  structure  This  basic  assumption,  combined 
with  Jacobi's  root  perturbation  formula,  leads  to  a 
fundamental  LAC  formula  for  predicting  algebrai¬ 
cally  the  root  shifts  produced  by  introducing  an  LAC 
structural  control  system  Specifically,  for  an 
undamped,  open- loop  structure  the  predicted  root 
shift  (dX^)p  is  given  by 


(dX^)p  '  2  ^ar  '’an  'rn  (II  1) 

where  the  coefficienl  matrix  C„,  is  a  matrix  of 
(damping)  gains,  and  d>ap.  <>,„  denote,  respectively, 
the  values  of  the  nth  mode  shape  at  actuator  station 
a  and  sensor  station  r 

Equation  (II- 1)  may  also  be  used  to  compute  the 
unknown  gains  C,,  if  the  dX„  are  considered  to  be 
desired  root  shifts  or,  equivalently,  desired  modal 
dampings.  While  an  exact  “inversion"  of  Eq  (II- 1) 
does  not  generally  exist,  weighted  least-squares 
type  solutions  can  be  devised  to  determine  the  actu¬ 
ator  control  gains  C,,  necessary  to  produce  the 
required  modal  damping  ratios.  This  determination 
of  the  gains  is  the  synthesis  of  LAC  systems. 

For  structures  which  already  have  some  damping 
or  control  systems  in  which  sensor,  actuator,  or  filter 
dynamics  can  either  be  ignored  or  are  already 
embedded  in  the  plant  dynamics,  the  root  perturba¬ 
tion  techniques  and  cost  function  minimization 
methods  above  can  similarly  be  used  to  synthesize 
low-authority  controls 

Robustness  of  LAC  systems.  When  sensors  and 
actuators  are  colocated  (i  c  .  a  =  r).  arc  comple¬ 
mentary.  and  only  rate  feedback  is  used,  formula 
(II- 1)  reduces  to 


which  shows  that  the  root  shifts  are  always  towards 
the  left  of  the  j-axis  if  all  the  gains  are  negative.  This 
robustness  result  is  obviously  based  on  the  assump¬ 
tion  that  both  sensors  and  actuators  have  infinite 
bandwidth,  and  also  that  the  structure  was  initially 
undamped.  Several  departures  from  this  idealiza¬ 
tion  occur  in  the  actual  practical  implementation  of 
the  LAC  systems  The  most  severe  of  these  results 
from  the  finiteness  of  the  actuators'  bandwidths 
More  precisely,  the  second-order  roll-off  introduced 
by  the  actuator  dynamics  will  always  destabilize  an 
undamped  structure.  However,  when  some  natural 
damping  is  present  in  the  structure,  or  when  a  pas- 


siv*.-  dnmpci  is  niouniod  in  parallel  with  thf  actua- 
toi,  additional  actu'i'  damping  can  bt  obtained  with 
out  destabilizing  the  structure 

High-Authoriry  Control  Design.  The  HAC  con¬ 
trol  design  procedure  can  be  based  on  any  number 
of  multivariable  design  methods,  e  g  .  LQG. 
Q- parameterization.  H*-optjmization.  etc  Increased 
penalties  in  the  LQG  cost  functional  are  placed  at 
those  frequencies  where  less  response  is  desired.  The 
concept  of  frequency-shaped  cost  functionals  was 
introduced  prior  to  ACOSS  (Gup.  1). 

The  frequency  shaping  methods  are  useful  in  sev¬ 
eral  areas  of  large  space  structures  control.  Three 
principal  applications  are  important:  (1)  robustness 
(spillover  avoidance).  (2)  disturbance  rejection,  and 
(3)  state  estimation 

Management  of  Spillover.  Spillover  in  closed- 
loop  control  of  space  structures  is  managed  by 
injecting  minimum  control  power  at  the  natural  fre¬ 
quencies  of  the  unmodeled  modes.  Procedures  for 
controlling  spillover  at  high-frequencies  are  usually 
discussed,  although  similar  techniques  are  applica¬ 
ble  for  other  regimes 

The  high-frequency  spillover  may  be  controlled 
by  modifying  the  state  or  the  control  weighting. 
Conversion  to  the  frequency  domain  gives  the  fol¬ 
lowing  performance  index: 


The  problem  of  robustness  (spillover  manage¬ 
ment  is  solved  by  making  Q  and  R  functions  of  fre¬ 
quency  Figure  11.2-36  depicts  the  modification  to 
the  nominal  LQG  controller.  Observe  that  fre¬ 
quency  shaping  adds  filters  whose  inputs  are  the 
innovations  outputs  of  the  state-estimator  in  the 
LQG  controller 


Fig.  II. 2-36  LQC  Control  With 

Frequency-Shaping  Fitters 

Summary.  The  application  of  frequency-shaping 
methods  to  large  space  structures  leads  to  a  linear 
controller  with  memory.  However,  additional  states 
are  needed  to  represent  frequency-dependent 
weights,  hence,  there  is  an  increase  in  the  controller 
order.  The  software  needed  for  these  controller 
designs 'S  similar  to  that  for  standard  LQG  problems 
Controller  Design  Using  Q-parametrization 
and  H'Optimization.  During  the  last  decade, 
mathematical  theories  of  servo  design  have  been 
based  mainly  on  quadratic  minimization  of  the 
Weiner-Hopf-Kalman  type,  usually  applied  to  state- 
space  models,  e  g..  LQG  controls.  However, 
despite  the  academic  success  of  these  methods, 
classical  frequency  response  techniques  relying  on 
“lead-lag  compensators"  to  reduce  sensitivity  have 
continued  to  dominate  industrial  servo  design  One 
reason  is  that  quadratic  design  tends  to  have  poor 
sensitivity.  On  the  other  hand,  the  frequency 
domain  description  has  proven  to  be  more  suitable 
to  characterize  uncenainties  which  arise  in  the  plant 
approximation  identification,  and  frequency 
domain  technique  usually  results  in  more  robust 
design,  e.g  .  frequency-shaped  LQG  can  be  viewed 
as  an  indirect  frequency-domain  design  approach 
Two  direct  multivariable  frequency  domain 
design  techniques  have  become  popular  in  recent 
years,  the  Q-parametrization  technique  and  the  H,- 
optimal  sensitivity. 


Q-Parametrization  Design.  Consider  the  linear 
unite  feedhack  systems  shown  In  Fig  II  2-37  where 
P(s)  IS  the  given  lii,».  me  invariant  plant.  C(s)  is 
the  linear  compensator.  U|  is  the  reference  input. 
Uo  and  d^,  are  respectively  the  plant-input  distur¬ 
bance  and  plant-output  disturbance,  and  yj  I's  the 
plant  output 


Fig.  11.2-37  The  Unity-feedback  Systems 


The  closed-loop  system  input-output  transfer 
function  is  given  by 

C  (t  PC)"’  PC  (1  I-PC)'’  -C  (I  PC)"' 

H  * 

[pC  (I  ♦  PC)’’  P  (I  CP)’’  (I  -.  PC)’’ 

(For  simplicity,  we  drop  the  argument  s  in  P(s). 
C(s).  etc.  in  this  section.) 

By  introducing  the  parameter  (transfer  function) 


0  =  C  (I  ♦  PC)’’ 
can  be  rewritten  as 

Q  -QP 
PQ  (I'PQIP 

No’e  that  the  closed-loop  input-output  transfer 
function,  for  the  given  plant  P.  is  completely  speci¬ 
fied  by  the  parameter  Q  in  a  very  simple  manner:  it 
involves  only  sums  and  products  of  P  and  Q 
In  a  typical  control  system  design  problem,  the 
two  most  important  closed-loop  transfer  functions 
are  and  H  is  the  transfer  function 

from  reference  input  u  j  to  output  y^  and  is  the 
transfer  function  from  plant-output  disturbance  d^, 


and  output  y_,  Ttiey  specify  respectivelv  the  servo- 
performance  and  regulator  performance  of  the 
feedback  system  S  The  two  transfer  functions  are 
given  by 

H  PQ  AND 


Therefore,  the  control  design  problem  reduces  to 
choosing  the  parameter  Q  so  that  the  closed-loop 
svstem  S  is  stable  and  that  „  and  H,  are 
"satisfactory"  After  the  parameter  Q  is  chosen,  the 
corresponding  compensator  C  can  be  obtained  by 
the  formula  (Cal  1.  Chap.  8] 

C  Q(i-PQ)  ’ 

Hence,  there  is  a  one-to-one  correspondence 
between  C  and  Q.  Consequently,  for  each  parame¬ 
ter  Q  chosen,  there  is  a  unique  compensator  C 
which  achieves  the  specified  Q 

The  selection  of  the  parameter  Q  in  the  design 
process  raises  several  questions:  What  arc  the  con¬ 
ditions  on  Q  so  that  the  resulting  compensator  C  is 
realizable  (e  g.,  proper)?  What  arc  the  class  of  all 
Q’s  which  result  in  a  stable  feedback  system?  How  is 
an  "optimal"  Q  chosen? 

Realizability.  If  the  plant  P  is  realizable,  then  the 
compensator  C  is  realizable  if  and  only  if  the  param¬ 
eter  Q  is  realizable.  Note  that  a  physical  plant  is 
always  realizable 

Global  Parametization.  If  the  open-loop  plant  P 
IS  stable,  then  the  closed-loop  system  S  is  stable  if 
and  only  if  Q  is  stable  since  sums  and  products  of 
stable  transfer  function  matrices  are  stable  Conse¬ 
quently.  the  class  of  all  stabilizing  compensators  is 
given  by 

|0  (I  PQ)  ’  0  IS  STABLE) 

and  the  class  of  all  achievable  stable  input-output 
transfer  matrix  and  the  class  of  all  achievable 
stable  disturbance-to-output  transfer  matrix  H 
are  given  respectively  by 

jPQ  I  Q  IS  STABLE  ) 

and 


(l  P0|  0  IS  STABLE) 


These  sets  give  global  parametrization  of  all  stabili- 
zaiing  compensators,  and  all  achievable  I '  O  charac¬ 
teristics  in  terms  of  a  stable  proper  transfer  matrix  Q 
In  other  words,  the  class  of  all  "feasible  "  designs  are 
parametrized  by  Q 

If  the  open- loop  plant  P  is  not  stable,  additional 
constraints  have  to  be  added  to  the  choice  of  Q.  in 
addition  to  stability  and  realiabiiity  of  Q.  For  exam¬ 
ple.  Q  must  contain  right  half  plane  zeroes  to  cancel 
the  unstable  poles  of  P  Currently,  there  are  three 
approaches  to  obtain  global  parametrization  of  a 
given  unstable  plant;  (i)  Factorization  representation 
theory  jDes.  11;  (ii)  Direct  approach  (Zam.  I); 
(iii)  Two-step  compensation  [Zam.  1). 

Optimality.  The  Q-parametrization  alone  does  not 
quantatively  address  the  issue  of  optimal  design. 
The  designer  selects  Q  from  the  class  of  "feasible" 
designs,  on  the  basis  of  the  desired  input-output 
response,  a  priori  knowledge  of  external  distur¬ 
bances.  bandwidth,  dynamic  range  and  uncertainty 
of  the  plant,  etc. 

Optimal  design  based  on  the  Q-parametrization 
and  fractional  representation  framework  has 
become  very  popular  in  the  research  community. 
The  Hoc-optimal  sensitivity  design  are  among  the 
results  available. 

Hoc. Optimal  Sensitivity  Design.  The  Hoc- 
optimal  sensitivity  design  is  an  extension  of  the 
Q-par;’  netrization  technique  to  include  a  quantita¬ 
tive  performance  measure  of  the  closed-loop  system 
and  achievable  optimality  based  on  the  perfor¬ 
mance  measure.  Roughly  speaking,  the  Hoc  design 
problem  is  the  following:  given  an  open-loop  plant 
Pis)  and  a  low-pass  weighting  function  W(s).  find 
the  compensator  C(s)  so  that  the  Hw-norm  of  the 
weighted  sensitivity  (1  +  PC)  ~  'W  is  minimized  sub¬ 
ject  to  the  stability  of  the  closed-loop  system 

Using  the  Q-parametrization  formulation,  the 
problem  is  equivalent  to  the  following:  find  a  Q  in 
Hoc  such  that  the  closed -loop  system  is  stable  and 
that  (I-PQ)W  is  minimized  Since  the  weighted 
sensitivity  function  is  affine  in  Q,  the  equivalent 
problem  is  easier  to  solve  than  the  original  problem. 
Solution  to  the  Hoc-optimal  Sensitivity  Prob¬ 
lem.  Based  on  the  fractional  representation 
(Coprime  factorization)  formulation,  several  solu¬ 
tions  have  been  proposed  and  algorithms  given. 


However,  all  the  proposed  algorithms  are  concep¬ 
tual  in  nature,  suitable  only  for  simple  text  book 
example.  More  effort  is  needed  toward  a  numeri¬ 
cally  robust  synthesis  procedure.  Relevant  research 
results  available  during  the  JOSE  period  of  perfor¬ 
mance  will  be  evaluated  and  used  as  appropriate. 
iI-2.2.7.4.2  Adaptive  Control  Techniques. 
Uncertainties  in  both  disturbance  spectra  and  sys¬ 
tem  dynamical  characteristics  will  limit  the  perfor¬ 
mance  obtainable  with  fixed  gain,  fixed  order  con¬ 
trols.  The  use  of  adaptive  type  control,  where 
disturbance  and  or  plant  dynamics  are  identified 
prior  to  or  during  control,  give  system  designers 
more  options  for  minimizing  the  risk  in  achieving 
performance  benchmarks.  For  the  case  of  SBL 
spacecraft,  where  performance  levels  are  extremely 
high,  it  is  absolutely  necessary  that  disturbance  and 
plant  models  be  equally  accurate.  Since  data  from 
ground  tests  do  not  usually  represent  the  flight  con¬ 
dition  accurately,  it  follows  that  an  on-line  proce¬ 
dure  for  identification  and  control  is  necessary. 


In  this  section  a  method  for  on-orbit  identification 
and  control  of  flexible  spacecraft  referred  to  as 
"adaptive  calibration",  is  described.  This  method  is 
being  developed  by  ISI  in  an  on-going  basic 
research  program  in  adaptive  control  supported  by 
AFOSR  Directorate  of  Aerospace  Science.  The 
basic  objective  of  this  research  program  is  to  estab¬ 
lish  the  theoretical  foundations  and  performance 
limitations  for  adaptive  control  applications  to  large 
space  structures.  An  important  element  of  the 
research  is  to  examine  implementation  concepts 
which  can  lead  to  appropriate  hardware  develop¬ 
ment.  A  summary  of  recent  results  is  in  [Kos  1] 
The  program  was  originally  formulated  in  late 
1982  in  response  to  the  increasing  concern  that  per¬ 
formance  robustness  of  Air  Force  LSS  type  system 
would  be  inadequate  to  meet  mission  objectives 
The  need  to  identify  modal  frequencies,  for  exam¬ 
ple.  in  high-performance  disturbance  rejection 
systems  has  been  shown  in  ACOSS  (1981)  and 
VCOSS  (1982).  The  deployment  of  high- 
performance  optical  or  RF  systems  may  require  on¬ 
line  identification  of  critical  medal  parameters  before 
full  control  authoritv  can  be  exercised  Parameter 
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scnsitiviiv.  manilfsifd  bs  porformanri’  degradation 
or  loss  o(  stability  (poor  lobustnesst  may  be  effec¬ 
tively  reduced  by  adaptive  feedback  mechaniza¬ 
tions  Reducing  the-  effects  of  on-board  disturbance 
sources  on  the  system  performance  (disturbance 
rejection)  is  particularly  important  for  planned  Air 
Force  missions  For  these  cases,  adaptive  control 
mechanizations  are  needed  to  produce  the  three-to- 
five  orders-of-magnitude  reductions  in  line-of-sight 
jitter  required  by  the  mission. 

AFOSR  Research  Program  Summary.  The 
research  was  originally  directed  toward  real-time 
adaptation  using  standard  estimator  and  controllers 
forms  in  the  basic  adaptive  control  structure,  as 
shown  in  Fig  11.2-38  Most  adaptive  control  algo¬ 
rithms  can  be  described  in  this  form  For  example, 
one  could  select  from  the  following  catalogs  of 
major  areas 


Model  Control  Design  Adaptation 

.ARMAX  Model  Reference  Gradient 

State-  Self-Tuning  Recursive  Least 

Squares 

Space  Pole-Placement  Recursive  Max 

Likelihood 
Extended  Kalman 
Filter 


Fig.  11.2-38  Adaptioe  Control  System 


The  schemes  also  differ  in  terms  of  update  rates. 
Typically  the  outer  control  loop  is  at  a  fast  rate, 
whereas  the  parameters  from  identification  are 
updated  more  slowly  Adaptive  schemes  are 


referred  to  as  "recursive  if  the  identification  rate  is  a 
fixed  multiple  r.f  the  controller  rate  If  identification 
is  used  w'hen  necessary  for  calibration  the  scheme  is 
referred  to  as  adaptive  calibration 

Although  a  great  deal  of  research  results  are  avail¬ 
able  about  adaptive  control  and  identification, 
unmodelcd  dynamics  and  broadband  disturbances 
will  significantly  upset  most  algorithms  Hence,  the 
lack  of  a  well- developed  robustness  theory  for  adap¬ 
tive  mechanizations  required  a  reexamination  of  the 
problem  at  a  more  fundamental  level.  i  e  .  develop¬ 
ment  of  model  and  disturbance  uncertainty  bounds 
for  which  adaptive  algorithms  would  exhibit  (stable) 
desired  performance  Toward  this  end  there  have 
been  two  major  accomplishments: 

(1)  Development  of  Theory  In  examining  the 
possible  use  of  recursive  adaptive  control  it  was  nec¬ 
essary  to  generate  new  theory  of  use  on  large  space 
structures  [Kos  3).  This  theory  accounts  for  the 
effect  of  unmodeled  dynamics  with  distributed 
parameter  systems,  such  as  flexible  space  struc¬ 
tures.  and  extends  current  adaptive  theory  in  sev¬ 
eral  directions 

In  the  first  place,  current  adaptive  theory  provides 
conditions  for  "global"  stability,  i.e..  bounded-input, 
bounded-output  stability  with  no  limitation  on  the 
size  (or  spectrum)  of  the  bounded-inputs  (e  g.,  dis¬ 
turbances  and  references).  Secondly,  the  theory  is 
limited  to  finite-dimensional  linear  systems  This  lat¬ 
ter  condition  cannot  be  satisfied  by  a  flexible  space 
structure,  which  is  a  distributed  parameter  system. 
Also,  the  disturbance  and  reference  inputs  effecting 
the  spacecraft  have  limited  magnitudes  and  spec- 
trums  and  these  limits  are  known,  although  not  pre¬ 
cisely  The  developed  theory  circumvents  those  dif¬ 
ficulties  by  providing  conditions  for  "local"  stability, 
i.e  .  limitations  in  input  size  and  spectrum  are 
accounted.  The  theory  also  allows  for  a  distributed 
system  as  well  as  providing  quantifiable  bounds  on 
permissible  model  error.  These  results  extend  the 
state-of-the-art  in  adaptive  theory  beyond  the  cur¬ 
rent  limits 

(2)  Adaptive  Calibration:  The  use  of  "slow" 
adaptive  control,  which  is  more  practical  than  recur¬ 
sive  adaptive  control  in  most  space  applications, 
necessitated  a  new  methodology  development 
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merging  kev  ideas  in  parameter  estimation,  system 
identification,  and  robust  conliol  design  The  term 
"slow"  means  that  there  is  sufficient  time  to  run 
batch  identification  before  the  control  system  is 
modified  The  methodology  developed  resolves  a 
long  standing  problem  with  adaptive  systems  of  this 
type,  namely,  the  means  to  provide  a  guaranteed 
level  of  performance  given  an  "identified"  model  of 
the  system  together  with  the  model  error  between 
the  system  and  the  identified  model.  In  fact,  the 
methodology  generates  performance  versus  model 
error  tables  {to  be  stored  in  the  computer)  from 
which  the  control  design  is  immediately  obtained. 
Moreover,  the  order  of  the  control  design  is  deter¬ 
mined  strictly  on  the  basis  of  model  error  and  per¬ 
formance  demand,  rather  than  trial  and  error  as  has 
been  suggested  in  the  past 

Application  of  Adaptive  Calibration.  The  basic 
problem  with  control  based  on  identified  models  is 
that  without  a  measure  of  model  error  it  is  very  easy 
to  dcstabili2e  the  system— particularly  when  the  goal 
is  high  performance — as  in  SBL  Adaptive  calibra¬ 
tion  is  an  approach  which  incorporates  a  measure  of 
model  error  with  robust  control  design  in  an  iterative 
way  so  that  identification  is  performed  only  where  it 
is  needed.  A  proposed  adaptive  calibration  system 
is  shown  in  Fig.  It. 2-39.  with  test  results,  using  the 
CSDL  #2  model,  shown  in  Fig  11.2-40.  The  adap¬ 
tive  calibration  procedure  involves  the  following 
steps: 

1.  The  model  M(s)  is  a  10-mode  model  which  has 
been  obtained  from  I/O  data 


PSD 


CSDL  MODEL  NO.  2 


Fig.  JI-2-40  Draper  Simulation  System 

2  Estimate  6(u;)  =  model  error  versus  frequency 
using  FFT.  This  is  dashed  curve  in  Fig.  II. 2-41. 

3.  Using  the  identified  model  M(s)  and  the  model 
error  6(w) .  synthesize  a  robust  control  (e.g. .  Sec¬ 
tion  2. 2. 7. 5). 

4  Calculate  —  stability  margin.  This  is  the 
dashed  curve  in  Fig. 11, 2-41.  Compare  to  model 
error  6  both  plotted  in  Fig  11.2-42.  If  acceptable 
go  to  Step  7  and  implement  controller.  Other¬ 
wise  go  to  Step  6. 

6.  Modify  filter  windows,  number  of  parameters 
(e  g.,  number  of  modes),  or  input  spectrum  and 
then  repeat  Step  1  to  obtain  new  ID  model. 
Fig.  II  2-42  shows  result  of  identification  after 
one  mode  is  added  in  the  frequency  domain 
region  where  the  test  fails. 


Fig.  11.2-39  Adaptive  Calibration  for  LSS 


Fig.  11.2-41  Comparison  of  Margin  and  Model 
Error 


Fig.  11.2-42  Data  Before  and  After  ID  Cycle  Plots 
Show  Where  Error  is  Margin 


7  Implement  controller 

This  strategy  has  been  successfully  implemented 
in  ground  test  hardware  Similar  techniques  arc 
being  used  for  adaptive  (hardware)  gun  turret  stabi¬ 
lization  for  an  Army  helicopter  weapons  platforms. 
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ABSTRACT 


A  method  of  averaging  is  developed  for  the  stability  analysis  of 
linear  differential  equations  with  small  time-varying  coefficients  which  do 
not  necessarily  possess  a  (global)  average.  The  technique  is  then  applied 
to  determine  the  stability  of  a  linear  equation  which  arises  in  the  study  of 
adaptive  systems  where  the  adaptive  parameters  are  slowly  varying.  The 
stability  conditions  are  stated  in  the  frequency-domain  which  shows  the 
relation  between  persistent  excitation  and  unmodeled  dynamics. 
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1.  INTRODUCTION 


For  a  large  class  of  adaptive  feedback  systems,  as  well  as  for  some 
output  error  identification  schemes,  a  stability  analysis  in  the  neighbor¬ 
hood  of  the  desired  behavior  leads  to  investigating  the  stability  of  the 
following  homogeneous  linear  system  of  differential  -  operator  equations 
(see  e.g.  [1]  -  [3]) 

e  =  -  c  u(t)  H(u(-)'e(  )).  Vt  e  (1.1) 


where  6(0)  e  R^ ,  c  is  a  positive  constant,  u("):  R^  ■*  R^  is  regulated  and 
bounded,  and  H  is  a  linear  time-invariant  operator  whose  transfer  function 
H(s)  is  proper,  rational,  and  stable,  i.e.,  all  poles  have  negative  real 
parts. 


Linearization  and  Local  Stability 


In  [2],  for  example,  system  (1.1)  is  obtained  as  a  result  of  lineari¬ 
zation  of  the  adaptive  system  in  the  neighborhood  of  a  "tuned”  system,  i.e., 

P 

a  system  where  the  adaptive  parameters  are  set  to  a  constant  value  6^  c  R 
and  whose  behavior  is  deemed  acceptable.  Hence,  in  (1.1),  e(t)  is  the 
vector  of  parameter  errors  between  the  parameter  estimate  at  time  t  and 
the  tuned  value  0^,  u(t)  is  the  regressor  vector  from  the  tuned  system 
(e.g.,  filtered  revisions  of  measured  signals),  and  the  scalar  e  is  the 
magnitude  of  the  adaptation  gain  which  essentially  controls  the  rate  of 
adaptation.  The  operator  H  depends  on  the  actual  system  being  controlled 
or  identified  and  also  on  the  tuned  parameter  setting  6,. 

It  is  shown  in  [?,3]  that  if  the  zero  solution  of  (1.1)  is  uniformly 
asymptotically  stable  (u.a.s),  then  the  adaptive  system  is  locally  stable. 


i.o.,  the  ndnptivc"  syst.oin  brlnvlot-  will  remain  in  a  neighborhood  of  the 
desired  behavior  provided  'he  initial  parameter  error  OfO)  and  the  effect  of 
external  disturbances  are  sufficiently  small.  Although  the  results  in  [?,3] 
were  arrived  at  using  input-output  properties  [16],  the  local  stability 
property  also  follows  from  the  results  on  "total"  stability  [i|]. 


Unmodeled  Dynamics  and  Slow  Adaptation 

In  the  ideal  case  there  are  a  sufficient  number  of  adaptive  parameters 
(the  number  p)  such  that  the  tuned  parameter  setting  results  in  H(s)  being 
strictly  positive  real  (SPR),  i.e..  Re  H(jii))  >  0,  Vw  c  R^ .  Under  these 
conditions,  we  have  the  following  results  (see  e.g  ,  [5]~[8],  [1]):  (1)  the 

zero  solution  of  (1.1)  is  stable,  i.e.,  e(t)  is  bounded  but  not  necessarily 
constant;  (2)  if,  in  addition,  u(t)  is  persistently  exciting,  then  the  zero 
solution  is  uniformly  asymptotically  stable  (u.a.s,),  thus,  e(t)  0 
exponentially  fast  as  t  •♦  “.  The  trouble  starts  when  there  are  an 

A 

insufficient  number  of  parameters  to  obtain  H(s)  e  SPR,  as  is  the  case  in 
adaptive  control  when  the  plant  has  unmodeled  dynamics  (see  e.g.  [2,  7], 
[12]). 

In  this  paper  we  will  examine  the  stability  of  (1.1)  when  e  is  small, 
u(t)  is  persistently  exciting,  and  H(s)  is  not  necessarily  SPR  but  only 
stable.  Reidle  and  Kokotovic  [9]  refer  to  this  case  as  "slow  adaptation" 
and  by  using  the  methods  of  averaging  described  by  Hale  [10],  they  show  that 
the  stability  of  the  zero  solution  of  (1.1)  is  critically  dependent  on  the 
spectrum  of  the  excitation  in  relation  to  the  frequency  response  H(ju)). 

With  the  same  assumptions,  Astrcm  [11]  uses  averaging  techniques  to  analyze 
the  interaction  between  unmodeled  dynamics  and  external  inputs  in  the 
counter-example  posed  by  Rohrs  et  al.  [12].  Both  these  analyses  require  the 
assumption  that  u(t)  is  almost  periodic.  In  this  case  Reidle  and 
Kokotovic  [9]  show  that  the  zero  solution  of  (1.1)  is  u.a.s.  If 


(1.2) 


-■ 


I  [a((»))a(u))*]  Re  H(Ju))>  0 


wli-'To  U  and  (.Ju),  (n  ■  Jil  ai'c,  respocti vely,  the  Fourier  exponents  and 
coefficients  of  u(t).  Condition  (1.2)  can  be  considered  as  a  post tl vi ty 
condi 1 1  on .  but  unlike  the  SPR  condition  Re  H(Ju))  is  not  required  to  be 
positive  at  al 1  frequencies. 

Averaging: _ Uses  and  Limi tations 

The  main  contribution  of  this  paper  is  to  extend  the  theory  of 
averaging  to  include  the  case  when  u(t)  does  not  have  a  (general Ized) 

Fourier  series  representation,  but  is  only  known  to  be  regulated  and 
bounded.  Thus,  u(t)  need  not  be  almost  periodic  nor  even  possess  a  (global) 
average  value.  We  also  state  stability  conditions  in  the  frequency-domain 
in  a  form  similar  to  (1.2).  Analagous  results  are  stated  for  the  discrete¬ 
time  system 

e(t+1)  -  0(t)  -  e  u(t)  H(u(«)'e(-)),  Vt  e  (1.3) 
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where  we  only  require  u(*)  e  and  H  to  be  1 inear-time-invariant  and 
stable.  Averaging  results  for  (1.3)  with  H  -  1  and  with  u(*)  not  almost 
periodic  can  also  be  found  in  [13];  and  this  suggests  the  possibility  of 
being  able  to  dispense  with  the  almost  periodicity  assumption  on  u(’)  and 
analyzing  (1.1)  with  a  non-SPR  operator  H. 

The  averaging  theory  developed  here,  as  well  as  averaging  theory  in 
general,  has  its  uses  and  limitations  for  adaptive  system.  In  the  first 
place,  the  theory  requires  slow  adaptation  which  can  be  counter-productive 
because  performance  can  be  below  par  for  the  long  period  of  time  it  takes 
for  the  parameters  to  readjust.  Secondly,  averaging  theory  is  a  form  of 
linearization,  thus,  the  (nonlinear)  adaptive  system  must  be  initialized  in 
a  (not  necessarily  small)  neighborhood  of  the  tuned  system.  On  the  positive 
side,  however,  we  do  obtain  frequency  domain  conditions  which  explain  the 
system  behavior  near  the  tuned  solutions.  In  this  sense,  we  can  consider 


thf  reyults  of  averaging  t  h''or-y  to  bo  necnssary  conditions  for  good 
performance  of  adaptive  systems. 

To  obtain  the  heralded  goal  of  frequency  domain  stability  conditions, 
it  may  be  inevitable  to  encounter  linearization.  Somewhat  less  intuitively 
appealing  results  can  be  obtained  without  resorting  to  direct  linearization 
or  averaging,  e.g.,  in  [?,3],  [1*)]  and  [15]  the  results  arise  from  a 
combination  of  small  gain  theory  and  perturbation  theory. 

Organization  of  Paper 

The  paper  is  organized  as  follows:  Section  P  develops  methods  of 
averaging  for  general  systems  described  by  linear  differential  equations  — 
both  homogenous  and  inhomogeneous  systems  are  considered.  The  reader  can 
rega’^d  this  section  independently  from  the  rest  of  the  paper,  because  the 
systems  of  linear  equations  considered  are  the  most  general  and  need  not 
arise  from  adaptive  systems.  In  Section  3  we  apply  the  general  results  of 
Section  2  to  (1.1)  and  obtain  frequency  domain  stability  conditions.  In 
Section  1)  we  analyze  the  effect  of  unmodeled  dynamics.  In  Section  5  we 
state  the  di screte- time  versions  of  the  results  obtained  in  Section  1,  and, 
as  in  Section  1,  these  results  are  of  general  interest. 


P.  METHOD  OF  AVERAGING 


?.1  LINEAR  HOMOGENEOUS  SYSTEMS 


In  this  section  we  will  consider  the  homogeneous  linear  system 


X  -  e  A(t)x  (2, 

Lemma  (2, 

Suppose  in  (2.1)  that  A(’)t  R^’'^  is  regulated  and  bounded. 

Vs,T  c  R^ ,  the  transition  matrix  4i(s+t,s)  of  (2.1)  is  given  by 


4i(s+T,s)  -  exp[cT  A^(s)]  +  R(s,eT) 


(2 


where 


A^(s) 


S+  T 

-  1  A(t)  dt 

s 


(2 


is  the  local  average  value  of  A(t)  on  the  interval  s  £  t  <  s  ♦  x,  and 

r(eT||Al|,)  (etIIaII^)*  exp(ex  |  [a]  |_)  (2 


See  Appendix  for  proof. 


Rorricjrks 


(1)  Assuming  that  A(t)  is  regulated  and  bound'='d  Is  sufficient  for  the 
existence  and  uniqueness  of  solutions  [17]. 

(2)  Observe  that  Lemma  (2.2)  is  valid  Vs,t  e  and  Vc  c  R.  In  the 
sequel  we  use  Lemma  (2.2)  only  for  the  case  when  e  >  0  and  ei  is  small. 


.  •  »  .  K  < 


The  stability  properties  of  (2.1)  can  be  established  by  application  of 
Lemma  (2.2)  as  stated  in  Theorem  (2.9)  below.  We  first  require 


Definition  [16] 


(2.6) 


The  function  ii(*):  C  R,  defined  by 


u(M)  =  lim  (|I  +  aMl  ■  1  )/a 
a+0  ^ 


(2.7) 


is  called  the  measure  of  the  matrix  M,  where  I*],  is  the  matrix  norm  on 
induced  by  the  vector  norm  j*|  on  c".  For  example,  if  I*]  is  the  Euclidean 
norm  then  u(M)  “  max  Xl(M+M*)/2].  For  any  norm  on  we  have  the  relation, 


-w(-M)  <  Re  X(M)  <  ij(M),  VM  e  C 


(2.8) 


I.  . 


I . n 


;:::d 


'  “•  *.  ■-  N,"*  V*  *«*  S 


Thrnrfm 


(?.9) 


Suppose  A't)  In  (2.1)  is  regulated  and  bounded  with  the  sequence  of 
local  average  values  {A.j,(kT),  Vk  e  Z^),  Then: 

(i)  If  0  T  '  0  and  a  >  0  such  that 


y[A.j,(kT)]  <  -  a,  Vk  e  2^  (2  10) 

then  5  n  >  0  such  that  VeT  e  (0,  n)  the  zero  solution  of  (2.1)  is  u.a.s. 

(ii)  If^T>0  and  a  >  0  such  that 

u[-A.j,(kT)]  <  -  a,  Vk  e  (2.11 ) 

then  3  n  >  0  such  that  VeT  e  (0,  n)  the  zero  solution  of  (2.1)  is  completely 
unstable. 


Remarks 

(1)  The  proof  (see  Appendix)  is  based  on  Lemma  (2.2)  and  the 
inequality  [16]: 

exp{-eTwr-A.j,(kT)])  <  i<|>(  (k+1  )T,kT)  -  R(kT,ET)| 

<  expIeT  p[A.j.(kT)]),  Vk  E  Z^  (2.12) 


7n 


wh^»re  and  R(’,-)  are  as  defined  in  Lemma  {?.?). 

(?)  Whenever  (?.10)  holds  we  have  |exp{  cTA^(  kT) )  |  <  1 ,  Vk  r.  which 
Insures  a  contraction  (for  small  cT)  on  the  Interval  s  £  t  £  s  T.  It  Is 
possible  to  weaken  condition  (?.10)  and  still  have  a  contraction  by  Just 
enforcing  lexpicT  A^(kT)}|  <  1  directly  as  is  done  by  Coppel  [18]. 

(3)  Note  that  Thereom  (2.9)  can  be  stated  in  terms  of  stronger 
conditions  on  u[A.j,(s)],  Vs  c  . 

Using  the  same  technique,  but  allowing  A{-)  (equivalently  A.j.("))  to 
possess  a  global  average,  we  obtain  the  following  sharper  result. 


Theorem  (2.13) 

Suppose  A(t)  in  (2.1)  is  regulated,  bounded,  and  has  a  (global) 

-  „nxn 

average  A  c  R  ,  i .e. , 


lim  A_(s)  »  A  (2.  Ill) 

T-><o  ‘ 


uniformly  Vs  e  with  Re  X(A)  *  0.  Under  these  conditions: 


(i)  If  ^  a  >  0  such  that 

Re  X(A)  <  -  a  (2.15) 

then  5  Eo  >  0  such  that  Ve  e  (0,  e„),  the  zero  solution  of  (2.1)  is  u.a.s. 


(ii)  If  ^  o  >  0  such  that 


max  R*  HA)  > 


« 


(?.16) 


then  >  0  such  that  Vi  c  (0,  e^),  the  zero  solution  of  (2.1)  is  unstable. 


Discussion 

The  results  in  Theorem  (2.9)  and  Theorem  (2.13)  generalize  some 
results  obtained  by  averaging  methods  such  as  those  described  by  Hale  [10], 
or  as  obtained  by  Coppel  [18]  using  the  notion  of  integral  smallness. 
Theorem  (2.13)  is  a  classical  result  of  averaging  theory,  except  that  as 
stated  allows  for  functions  which  are  not  necessarily  almost  periodic.  The 
class  of  functions  allowed  in  theorem  (2.13)  ■  regulated,  bounded,  with  an 
average  -  is  not  precisely  characterized.  Obviously  it  includes  the  class 
of  asymptotically  almost  periodic  functions  of  the  form  [19] 


A(t)  =  Ap(t)  +  A,(t) 


(2.17) 


where  ^^(t)  is  almost  periodic  and  A^(*)  e  . 

Theorem  (2.9)  considers  a  larger  class  of  functions  —  those  without 
an  average  —  at  the  expense  of  a  weaker  result:  the  stability  — 
instability  boundary  is  not  as  sharp  as  in  theorem  (2.13). 

An  example  of  a  function  which  satisfies  the  conditions  of  theorem 
(2.9),  but  not  of  (2.13)  is: 


A(t) 


Ao  * 


(I//2)  A^  (sin  log  t  +  cos  log  t) 


(2.18) 


where  A^  -  A^  >  0,  1«1,  2  such  that  A,  -  A^  >  0.  This  function  does  not 
have  a  global  average,  as  can  be  seen  from 


S"! 

^  1  A(t)dl  -  A„  ♦  A^rsin  ]og  (s^T)  *  |  sin  log  (1  -*  i)]  (?.19) 
s 

However,  it  satisfies  the  conditions  of  Theorem  (2.9)  because  from  (2.18), 


s+T 

Y  \  A(t)  dt  >  Ao  -  A^  >  0,  Vs  c  R^,  VT  >  0.  (2.20). 

s 


From  the  proof  of  Theorem  (2.9)  we  can  extract  a  value  for  n  and  also 
state  bounds  on  the  exponential  rates  of  growth  or  decay  of  the  transition 
matrix  (JiCt.  t)  for  all  t  >  t.  Specifically,  we  have: 


Corollary  (2.21) 

If  A(t)  is  regulated  and  bounded  with  |1a(*)||„^  m,  then: 

(i)  Whenever  (2.10)  holds  for  some  T  >  0,  the  zero  solution  of  (2.1) 


is  u.a.s. 

VeT  t  (0,  n)  where: 

(a) 

n  >  0  satisfies  exp  (-no)  +  r(nm)  -  1 

(2.22) 

(b) 

Vt,T  E  with  t  >  T, 

(2.23) 

1((i(t,T)|  <  K  exp(-e(t-T)B) 

where 

K  -  exp(cT(m+B))  >  1 

exp(-c  TB)  •  exp(-cTa)  +  rCeTm)  <  1 
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(ii)  Whonpvor  (2.11)  holds  for  somo  T  ■>  0,  the  zero  solution  of  (2.1) 
is  unstable  VcT  c  (0,n)  where: 


where 


(a)  n  >  0  satisfies  exp  (not)  -  r(nin)  -  1 


(b)  Vt,T  c  with  t  >  T, 


(4i(t,-t)|  2  exp(e(t  t)B) 


K  -  exp(-cT(m+B))  <  1 


exp(eT6)  -  exp(cTa)  -  r{eTm)  >  1 


(2.2i)) 


(2.25) 


•  .  *  “  •  ^  j 

7^  CP 

'•  k  w'“f  4^ 

•  I 

►  • 


2.2  LINEAR  INHOMOGENEOUS  SYSTEMS 


In  this  section  we  extend  the  results  of  Section  2  to  the  inhomo¬ 


geneous  system 


X  -  e(p(t)  +  A(t)x) 


Theorem 


(2.26) 


(2.27) 


i^VkV-k.--. 


■ 

-  LV-  .V 

.V 

.V  -V*  .‘-V 
-  •  .  •  ^  ^ 
k  *  ■  •  '  J.  •  K 


Suppose  in  (2.26)  that  A(*):  R"*"  and  p(*):  R^  r"  are  regulated 

and  bounded.  Let  A^(s),  Vs.t  e  R^  denote  the  local  average  value  of  A(t)  as 
defined  by  (2.ii)  and  let  p^(s)  denote  the  local  average  value  of  p(t) 
defined  by 


*■  »  -v'  V  • 


(s)  -  -  J  p(t )  dt. . 


(?.?8) 


Under  these  conditions,  Vs,t  e  , 

|x(s+T)-(i>(s+T.s)''x(s)  +  CTp^(s)]|  <  Q  ( c  T ,  |  |  A  |  |  _^ )  |  |  p  |  |  ^  (2,29) 


where 


q(cT.||A||^)  -  (£t)"| lA| |^exp{ET|  |a! 1^) 


(2.30) 


Combining  this  result  with  Theorem  (2.9)  and  Corollary  (2.21)  gives 


Corollary  (2.31) 

Under  the  conditions  of  Theorem  (2.27),  if  ^  a  >  0  and  T  >  0  such  that 


u[A.j.(kT)]  <  -  a,  Vk  E  (2.32) 

then  VeT  e  (0,  n) , 

I  1*1  L  1  'll  'IjHpI  L  (2.33) 

keZ 

♦ 


with 


q,  -  eT  exp(ET(m-6))/(1  €Xp(-eT8)) 


(2.34) 


q,  -  eT  exp(ETra)[1  +  cTm  exp(eTm)/(1  exp(-ETe)] 


whero  m  -  I|a1|^  and  whore  ti  >  0  and  B  >  0  satisfy  (2.??)  and  (?.?3). 
respectively,  l.e., 


exp(-noi)  +  (nni)’exp(nni)  •  1 

expC  cTB)  “  exp(-cTa)  +  (eTni)^exp(cTml 


(2.35) 


Remarks 


Observe  that  as  eT  -*■  0,  qi  -»■  1/6  and  Qj/eT  1.  Hence,  as  eT  ■*  0  we 
see  that  ||x||^  is  overbounded  by  the  largest  value  of  the  sequence  of  local 
averages  {|p^(kT)(,  k  e  Z^}  and  not  |(p|(„.  For  example,  if  p(t)  -  sin 
then  ||p||  =  1  whereas  |p_(kT)l  -  0,  Vk  e  Z  . 


.  FREQUENCY-DOMAIN  STABILITY  CONDITIONS 


In  this  section  we  apply  the  results  of  Section  2  to  the  homogeneous 
linear  system  (1.1),  i.e., 


0  -  -e  u(t)  H(u(-)'e(  •)).  0(0)  e 


(3.1) 


where  H  is  a  linear  time  invariant  operator  with  transfer  function  H(s). 
We  first  show  that  for  sufficiently  small  e  >  0,  the  stability  analysis  of 
(3.1)  can  be  determined  from  the  stability  of  an  "averaged"  system 


0  =  e  avg(u(t)(Hu)’ (t)}0 


where  avg  {•}  has  yet  to  be  precisely  defined.  Using  this  result  we  then 
establish  stability  conditions  in  the  frequency-domain  involving  the  Fourier 
transform  H{ju))  and  the  "spectral"  content  of  u(t),  where  this  notion  has 
also  to  be  defined.  Finally,  we  show  that  the  appropriately  defined 
spectral  content  of  u(t)  necessarily  requires  that  u(t)  have  a  persistency 
of  excitation  property,  and  that  the  dominant  excitation  should  be  at  those 
frequencies  where  Re  H(ju)  >  0. 

The  first  step  in  the  analysis  is  to  transform  (3.1)  into  a  form 
suitable  for  application  of  Theorem  (2.9). 


Lemm.i 


Assume  In  (3.1)  that  u(t)  Is  regulated  and  bounded,  and  H  is  a 
causa]  linear  time-invariant  operator  where  transfer  function  H(s)  is 
proper,  rational,  and  stable.  Under  these  conditions,  6  Co  ^  0  such  that 
Ve  e  (0,e^),  (3.1)  is  equivalent  to. 


e  -  -  E  rB(t)  +  e  u(t)v(t)’]e 


(3.3a) 


where 


B(t)  =  u(t)(Hu(-)’ ){t) 


(3.3b) 


and  where  the  function  v(*)  :  R^  is  regulated  and  bounded. 


Remarks 


(1)  Since  v(t)  is  bounded,  it  follows  that  we  can  approximate  (3.3) 
to  first  order  in  e  by  6  »  -  e  B(t)9  which  is  the  form  required  by  Theorem 
(2.9).  Precise  conditions  on  this  approximation  are  stated  in  Theorem 
(3.11)  below. 

(2)  Estimates  for  e^  and  (|v|(^  are  obtained  from  the  proof  of  Lemma 
(3.2).  Specifically,  let  H(s)  -  d  +  c’(sI-A)  ^b  where,  with  no  loss  in 
generality,  (A,b,c,d)  is  minimal,  (c|  «  1,  and  Re  X(A)  <  0.  Then: 

e,  -  a(l  -  te)/[Kllul|.(d||u|l„  *  i)l] 


l|v||.  <  (1  -  to/U/c 


(3.5) 


whc-re  K  1  and  a  >  0  are  found  from 

|exp(tA)|  <  K  exp  (-ta),  Vt  t 


with 


1 

£  -  £,[1  4  (1  ♦  d||ullyi„)' 

K  •  IlLolL 


where  L^Ct)  satisfies, 

L,  -  AL„  +  bu(t)’ ,  L„(0)  =  0  e  r"^^ 


(3.6) 


(3.7) 

(3.8) 

(3.9) 


We  now  use  Lemma  (3.2)  to  establish  frequency  domain  stability 
conditions  for  (3.1).  This  requires  that  u(t)  be  restricted  to  those 
functions  which  have  a  Fourier  series  representation  on  any  finite  interval. 
A  known  class  of  such  functions  is  defined  as  follows. 


Definition  [20] 

A  function  f('):  R  r"  is  a  c"  function  if  it  is  regulated,  bounded 
+  0 

and  ■]  a  constant  6  >  0  such  that  any  two  points  t,,  t^  e  R^  where  f(*)  is 
discontinuous  are  separated  by  at  least  an  interval  6,  i.e.,  |t,  ~  ^  6. 


Frequency-domain  stability  conditions  for  (3*1)  can  now  be  stated. 


Thf“Or-pm 


(3.11) 


Assumo  in  (3.1)  that: 

(Al)  H  is  linear  with  stable  proper  rational  transfer 

A 

function  H(s)  and  impulse  response  h(t).  Thus,  a  >  0  and  b  >  0  such  that 
|h(t)  -  h(0)6(t)|  <  a  exp(-bt),  Vt  e  (3.12) 

(A2)  u(*)  c  with  piece-wise  Fourier  series  representation, 

g 

jo)  t 

u(t)  -  I  ^  ,  Vt  e  Cl<T,(k+1)T],  Vk  e  Z_^  (3.13) 

meZ 

for  any  T  >  6  where  oj  «  2irm/T. 

—  m 

Define  the  matrix  sequence  in  R*^*^  by 

R  (k)  -  I  Vk  e  Z^  (3.1^) 

mcZ 

Under  these  conditions; 

(i)  If  ?  T  >  6  and  a  >  0  such  that 

v[-R.j.(k)]  <  -  (a  *  2(a/b*)nullyT),  Vk  e  (3.15) 

then  "B  Ej  >  0  such  that  Ve  e  (0,  e„)  the  zero  solution  of  (3.1)  is  u.a.s. 

(ii) IfBT>6  and  a  >  0  such  that 


ufR^CkT)]  <  -  ((,  ^  P(a/b  )||u||VT).  Vk. 


(3.16) 


then  c  >  0  such  that  Vc  c  (0,  £„)  the  zero  solution  of  (3.1)  is  unstable. 


Remarks 


(1)  The  existence  of  the  piece-wise  Fourier  series  representation 

(3.13)  for  u(t)  is  guaranteed  by  u(*)  e  [17].  The  Fourier  coefficients 

0 

a,  (w  )  are  the  coefficients  of  the  T-periodic  function 
km 


jw  t 

“k^^^  ‘  ^  ,  Vt  e  R_^ 

meZ 


(3.17) 


which  is  equal  to  u(t)  for  kT  _<  t  _<  (k+1)T  and,  in  general,  not  equal  to 
u(t)  on  any  other  interval.  Thus,  U|^(t)  is  just  u(t),  Vt  t  [kT,  (k+1)T], 
repeated  with  period  T.  Observe  that  the  spectrum  of  U|^(t)  is  what 
determines  the  stability-instability  boundary  and  not  the  spectrum  of  u(t). 
These  will  merge  only  when  u(*)  has  a  (global)  average  as  assumed  in  Theorem 
(3.23)  below. 

(2)  The  matrix  R.j.{k)  can  be  equivalently  expressed  as  the  local 
average  value  of  Uj^(t) (HUj^)^  (t),  l.e.. 


• .  • . 


v>, 

•.■.V’.V.' 


(k+1)T 

R  (k)  -  4  \  u  (t)(Hu  );(t)dt 

kT 


(3.18) 


»  • 


where  (Hu,  )  (t)  is  the  "steady-state"  part  of  (Hu,  ){t),  i.e., 
k  “  k 


(Huj^).(t)  -  I  H(Ju.^)a^(u,^)e  ,  Vt  t  R^ 


(3.19) 


.  •  ’ . ' .  ^ 
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(3)  If  we  use  the  mermure  ij(K)  -  max  then  (3.  IS)  and 

(3.10)  become,  respectively, 

X[Q^(l<)]  >  a  +  ?(a/b")|  |u|  I  VT,  Vk  e  (3.20) 

and 

A[Q.j.(k)]_<  -(a  +  2(a/b^)||u||VT),  Vk  e  (3.21) 

where  Q^(k)  is  the  Hermitian  part  of  R.j.(k),  i.e., 

Q.j.(k)  •=  I  Re[aj^(ai^)aj^(a>^)*]  Re[H(ja)j^)]  (3.22) 

meZ^ 

(5)  The  "initial  conditions"  at  t  -  kT  contribute  to  the  term 
2(a/b^) I [u] I VT  in  (3. 15)-(3.l6)  or  (3.20)-(3. 21) .  Hence,  the  average 
energy  in  U|^(t ) (HUj^)_^(t )  must  dominate  long  enough  (T  sufficiently  large)  to 
overcome  these  (possibly)  negative  effects. 

As  before,  if  u(t)  is  further  restricted  such  that  R^(0  has  a  global 
average,  then  we  can  sharpen  the  stability-instability  boundary.  For 
example,  if  u(t)  is  almoso  periodic  then  a  Fourier  series  representation 
exists  Vt  E  R^and  R.j,(’)  has  an  average  [10].  The  stability  conditions  for 
this  case  are  stated  as  follows. 


Thf'Oi’em 


where  £2  e  R  are  the  distinct  Fourier  exponents  and  {a(u)),  w  e  P)  are  the 
Fourier  coefficients.  Define  the  matrix  R  e  by 


R  =  5]  a((D)  a(to)*  H(-ja)) 

U)  E  D 


(3.25) 


If  Re  A(R)  *  0  then  ^  e^  >  0  such  that  Ve  e  (0,  e^)  the  zero  solution  of 
(3.1)  is: 

(i)  u.a.s  if  Re  X(R)  >  0  (3.26) 

(ii)  unstable  if  max  Re  X(R)  <  0  (3.27) 


Discussion 

The  proof  of  Theorem  (3.23)  is  entirely  analogous  to  that  of  Theorem 
(2.13).  Theorem  (3.23)  is  the  result  obtained  in  [9]  when  u(t)  is  almost 
periodic.  Theorem  (3.11)  is  a  generalization  in  that  u(*)  e  . 

Observe  that  the  stability-instability  boundary  determined  by  (3-20)  - 
(3.22)  exists  if  and  only  if 


AtQ.j.(k)]  -  0,  Vk  E  (3.28) 
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By  (3.??),  this  will  hc<l  d  if  and  only  if  for  somo  finite  integer  q  _>  P  f 

rank  rcX|^(0),  -  p ,  Vk  e  .  (3.29) 

Hence,  Theorem  (3.11)  implicitly  restricts  li(‘)  to  those  functions 

whose  (time-varying)  Fourier  coefficients  satisfy  the  rank  condition  above. 
This  class  of  functions,  however,  are  precisely  those  which  can  be  cate¬ 
gorized  as  persistently  exciting: 

Definition  [1]  (3-30) 

A  function  f(*):  — >  R*^  is  persistently  exciting  over  an  interval 

if  it  is  regulated,  bounded,  and  -3  constants  h  >  0  and  B  >  0  such  that 

s+h 

min  X(^  1  f(t)f(t)’dt)  >  B.  Vs  c  R^  (3-36) 

s 

Denote  such  functions  by  u(’)  e  PE^(h,8). 

Hence,  we  immediately  see  that  if  u(t)  in  (3.1)  is  in  PE'^(h,B)  and 
then  (3.29)  will  hold  for  VT  h  >  6.  It  is  important  to  emphasize, 
however,  that  even  if  u(t)  is  PE,  u.a.s.  of  the  zero  solution  of  (3.1)  is 
guaranteed  if  (3.20)  holds.  The  implication  then  is  that  u(e)  must  have  a 
dominant  spectrum  at  those  frequencies  where  Re[H(jio^)]  >  0.  Thus,  we  can 
view  (3.20)  as  a  generalized  positivity  condition  on  the  operator  H. 

Condition  (3.20)  is  significantly  weaker  than  the  usual  positivity 
conditions  on  H.  For  example,  a  strictly  proper  transfer  function  H(s)  is 
strictly  positive  real  (SPR)  if  it  is  exponentially  stable  and  3  constant 
p  >  0  such  that  [16]: 
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R<  rH(j,.o]  >  p1h(j...)|  . 


(3.3.'") 


This  condition  must  hoi d  at  every  frequency,  whereas  (3.20)  requires 
Re  rH(j(»i^)]  >  0  at  those  di screte  frequencies  in  where  the  magnitude  of 
the  input  spectrum  is  large.  Conversely,  at  those  frequencies  in  R^  where 
Re[H{ja)^)]  <  0,  the  magnitude  of  the  input  spectrum  should  be  small.  Since 
(3.25)  will  fail  if  Re  H(i(o)  <  0,  Vw  c  R^  ,  it  follows  that  Re  HCjtv)  >  0  at 
some  frequencies,  hence,  the  motivation  to  refer  to  (3.20)  as  a  positivity 
condition. 

Although  condition  (3.20)  is  weaker  than  condition  (3-3?).  we  do  pay 
the  Piper.  Suppose  H(s)  is  SPR  and  (3.32)  holds.  If  u(t)  is  persistently 
exciting  than  Theorem  (3.11)  states  that  the  zero  solution  of  (3.1)  is 
u.a.s.  for  sufficiently  small  e  >  0.  However,  from  other  arguments  (see 
e.g.  [1])  we  know  that  under  these  same  conditions  the  zero  solution  of 
(3-1)  is  u.a.s.  for  all  c  >  0.  Thus,  Theorem  (3.11)  is  conservative  in  this 
case.  However,  when  H(s)  is  not  SPR  but  (3.32)  holds  at  some  frequencies. 
Theorem  (3.11)  is  now  applicable  whereas  the  results  in  [1]  do  not  apply. 

In  fact  in  this  latter  case  when  e  gets  too  large  the  zero  solution  of  (31) 
can  be  unstable,  even  if  (3.20)  holds.  For  example,  if  in  (3.1)  u(t)  - 
sin(0.35t)  and  H(s)  -  1/(s*  +  2s  +  2)  then  condition  (3.20)  is  satisfied. 

The  simulations  in  Figure  3.1  with  6(0)  «  1  show  that  the  zero  solution  is 
u.a.s.  for  e  *  If  but  is  completely  unstable  for  e  -  8. 


sTaV,!^  Q  ^  -  y') 
1 1  =  S') 
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H.  EFFKCT  OF  UNMODELED  DYNAMICS 


In  this  section  we  will  consider  the  system 


0  -  -e  u(t)  H(u(-)’e(-)),  6(0)  E  rP 


('4.1) 


where  H  as  before  is  linear  with  stable  transfer  function  H(s).  In  addition 
we  assume  that  H(s)  has  the  decomposition 


H(s)  -  H„(s)  A(s) 


(".2) 


where  H^Cs)  is  SPR,  i.e.,  p  >  0  such  that 


Re  H„(ju))  >  p|H„(ja))|^  Vu)  e  R^ 


(".3) 


and  where  A(s)  represents  unraodeled  dynamics  such  that 


|A(Ju)|  <  6j^((d),  Voj  e  R^ 


('l.'O 


We  also  assume  that  u(t)  satisfies  the  conditions  of  Theorem  (3.11)  in  that 
9  T  >  0  such  that  u(t)  has  the  piece-wise  Fourier  series  representation 


Ju)  t 
fn 


kT  <  t  <  (k+1)T 


u(t)  "  I  o  (u)  )e 


(".5) 


w>io(’c  »  pTim/T  and  k  i  .  We  also  decompose 


-  a  -(u  )  +  B  _(lD  ), 
kT  m  kT  m  kT  m 


Vm,k  c  Z 


(i<.6) 


where  a,  _(u)  )  is  due  to  predetermined  inputs  and  6. -.(id  )  is  due  to  disturb- 
kT  m  kT  m 

ances  bounded  by 


(«.7) 


Hence,  the  functions  oi  -*  and  m  -►  represent,  respectively, 

bounds  on  the  effect  of  unmodeled  dynamics  in  H{s)  and  unknown  elements  of 
u(t)  as  a  function  of  frequency.  Combining  the  above  assumptions  with 
Theorem  (3.11)  and  using  (3.25)  gives; 


w 


t'.’ 


Lemma 


(i1.8) 


The  zero  solution  of  (i(.1)  is  u.a.s.  if  e  >  0  is  sufficiently  small 
and  if  3  a  >  0  and  T  >  0  such  that  Vk  e  , 


Ml  PlH,(Ju.„)p  Re[X».j(w^)])  > 


q  ♦  a  ♦  2(a/b*)||u||,/T 


('1. 9) 


I 


vS-:-:-: 


where 


kT '  m ' 


(‘t.10) 


X°  (u)  )  -  o._,(a)  )  a  _((!)  )* 
kT  m  kT  m  kT  m 


(i(.11) 


Discussion 


Condition  (^.6)  shows  that  the  dominant  excitation  must  act  in  the 
frequency  range  where  H(s)  »  H^Cs)  c  SPR.  Moreover,  there  must  be  enough 
excitation  and  positivity  in  Re  H„(ju))  in  this  range  to  overcome  initial 
conditions  (the  1/T  term)  and  the  effect  of  unmodeled  dynamics  and  unknown 
disturbances  (the  term).  Typically,  the  disturbances  and  unmodeled 
effects  occur  at  high  frequencies  and  the  known  efforts  in  Hp(s)  and 


at  low  frequencies.  For  example,  if  there  is  a  frequency  u>^  such  that 


a°.j.(a)),  Ho(j(o)  small  for  w  >  > 


6„(u)),  6  (u)  small  for  at  <  <  w 
H  ’  u  c 


then  condition  (iJ.8)  holds  if  Vk  e  Z^, 


(*<.12) 


M  I  p|H„(Ja.„)I*  Re[X«.^(a,^)])  >  >  I  6j^(u)^)6» (u,^)  (i<.13) 


Observe  that  robustness  conditions  (4.8)  or  (4.13)  are  dependent  on  the  in¬ 
put  signal  spectrum  as  as  well  as  the  unmodeled  dynamics.  In  non-adaptlve 
linear  systems  the  robustness  conditions  only  involve  system  dynamics. 


5.  DISCRETE -TIME  SYSTEMS 


In  this  section  we  state  the  di screte- time  versions  of  Lemma  (2.2), 
Theorem  (5.2)  and  Theorem  (2.13)  for  the  linear  homogeneous  difference 
equation 


x(t+1)  -  (I  +  c  A(t))  x(t),  t  E 


(5.1) 


The  results  are  identical  to  the  continuous- time  results,  and  the  same 
comments  apply  mutati s  mutandis.  The  proof  of  the  following  Lemma  is  in 
appendix. 


Lemma  (5.2) 

Suppose  in  (5.1)  that  A(-)  e  J,"’'".  Then  Vs,t  e  ,  the  transition 
matrix  (|i{s+t,s)  of  (5.1)  is  given  by 


4)(s+t,s)  »  I  +  et  A^(s)  R(s,et) 


(5.3) 


where 


A^(s) 


2 


S-f  T-1 
I  A(t) 


t-s 


(5.^) 


is  the  local  average  value  of  A(t)  on  the  interval  s  £  t  £  s  ♦  t,  and 


the 


I|R(-.  ”)IL  < 


-  n  *  ‘  >I|A|L,).  V.  >  0 


^(”I|A1L)^  Vc  <  1/1|A| 


(5.5) 


Remarks 


(1)  A  similar  result  can  be  found  in  [13]. 

(2)  The  conditions  on  A(»)  :  -►  are  weaker  than  those  imposed 

in  the  conti nuous-time  version  in  Lemma  (2.2).  In  the  discrete-time  case  we 
only  need  A{*)  e  whereas  in  conti  nuous-time  A(»)  is  regulated  and 

bounded. 


The  following  stability  result  follows  immediately  from  Lemma  (5.2). 


Theorem 


(5.6) 


Suppose  in  (5.1)  that  A(*)  e  with  the  sequence  of  local  average 

values  {A^(kT),  k  c  Z^).  Then: 


(i)  If  ^  T  e  Z  and  a  >  0  such  that 


li[A^(kT)]  <  -  o.  Vk  e  Z^ 


(5.7) 


then  •]  n  >  0  such  that  VeT  e  (0,  n)  the  zero  solution  of  (5.1)  is  u.a.s. 


(ii)  If^T  eZ^  and  a  >  0  such  that 


p[-A^(kT)]  <  -a.  Vk  e  Z^ 


(5.8) 


.■  '.v.v 
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then  -3  n  >  0  suoh  that  ViT  t  (0,  n)  the  zero  solution  of  (5.1)  Is  unstable 

If  we  let  A(*)  have  a  (global)  average  than  by  applying  the  same 
argument  in  the  proof  of  Theorem  (2.13)  we  obtain  the  following  analogous 
result. 

Theorem  (5.9) 

Suppose  in  (5.1)  that  A(*)  e  with  (global)  average  A  e 

i  .e. , 

limA_.(s)-A  (5.10) 

uniformly  Vs  £  where  Re  X(A)  •  0.  Under  these  conditions: 

(i )  If  5  a  >  0  such  that 

Re  X(A)  <  -  a  (5.11) 

then  5  e,  >  0  such  that  Ve  e  (0,  e^),  the  zero  solution  of  (5.1)  is  u.a.s. 

(ii)  lf^a>0  such  that 

max  Re  A(A)  2  ®  (5.12) 

then  ^  Eo  >  0  such  that  Ve  e  (0,  e,),  the  zero  solution  of  (5.1)  is 


unstable . 


APPENDIX 


A.  Proof  of  Lemma  (2.2) 


Using  the  Peano-Baker  series  representation  for  the  transition  matrix 
of  (2.1)  gi ves: 

s+  T 

(|)(s+T,s)  =  I  +  e  1  A(t)clt 
s 

s+ T 

+  I  1  A(t,)  j  A(t,) 

k«2  s  s 


k-1 


A(tj^)(jt,  .  .  .  dt^  (A.1) 


Using  definitions  (2.3)  -  (2.^1)  for  R(s,eT)  and  A^(s),  respectively, 
together  with  the  series  expansion  for  exp(eT  A^(s))  results  in, 


R(s,ct)  "  I  [-(et  a  (s))*^/k! 
k-2 


S+  T 


t. 


k-1 


+  £*'  1  A(t,)  /  A(tj)  .  .  .  /  A(t  )dt,, 


.  dt|^]  (A. 2) 


<  ?I  (et 
h-2 


j'^/k!  ,  Vs  E  R^ 


(A. 3) 


exp(ETllAll  ) 


(A. it) 


since  ||A^(*)|1„  ■  l|A(*)||,.  This  proves  (2.5). 


B .  FrooT  of  Thforom  (?. 9)  and  Cot'Oilfit-y  (?.?1) 

The  following  Inequality  holds  VM  c  and  Vt  c  [16]: 

exp[-tu(-M)]  <  |exp(tM)|  £  exp  [tu(M)].  (B.1) 

Combining  this  inequality  with  (2.3)  gives. 

exp[-ETii(-A_^(s)]  _<  l4i(s+T,s)  R(s,ct)| 

<  exprcTM(A^(s) )] ,  Vs,T  c  R^  (B.2) 

which  implies, 

j(ti(s+T.s)'  <  exp[eTu(A^(s))]  +  r(CTm)  (B.3) 

|4i(s+T,s)|  2  sxp[-ci:u(A^(s))]  r(ETm)  (B.il) 

where  we  have  used  (2.5)  with  ||Aj|^  »  m. 

We  first  prove  part  (i)  by  using  condition  (2.10)  and  inequality  (B.3) 
with  T  -  T  and  s  -  kT.  This  gives, 

|(K(k+1)T,  kT)|  <  exp(-ETa)  +  r(ETm),  Vk  e  (B.5) 

We  now  need  to  show  that  e(kT)  -►  e((k+1)T)  Is  a  contraction  mapping,  i.e., 
the  right  hand  side  of  (B-5)  la  less  than  one  for  sufficiently  small  cT, 
I.e.,  9  n  >  0  and  B  >  0  such  that  V  eT  e  (0,  n). 


exp(-iTu)  •*  r(.Tm)  -  exp{-iTB) 


(B.6) 


e  xp  ( -  nu )  *  r  ( nm )  -  1 


From  the  definition  of  r(-)  in  (?.5)  it  is  obvious  that  (B.6)  holds. 

Observe  that  the  expressions  in  (B.6)  appear  also  in  conditions  {2.22), 
{2.23)  of  Corollary  {2.2^  ) . 

For  any  t,  s  e  with  t  2  s,  there  exists  an  integer  (<2  0  such  that 
s+kT<t<s+  (k+1)T.  Thus, 


|<()(t.s)|  -  |(J>(t,s+kT)({,(s+kT,s+(k-1)T) 


i}>(s+T,  s)| 


2  I t ( t ,s+kT) j exp  (-ckTB), 

<  |(f.(t,s+kT)|exp  (-e(t-s-T)e), 

<  exp(cT(m+6) )exp  (-c(t-s)6) 


by  (B.6) 
by  kT  >  t-s-T 

(B.7) 


by  the  inequality  [16]: 


t 

|(J)(t,s+kT)|  <  exp(  1  u[EA(T)]dT) 

s+kT 


<  exp(Em(t-s-kT) ) , 


<  exp ( emT ) 


(B.8) 


by  ii[A(T)]  <  |A(t)|  2  m 


(B.9) 


by  t-s-kT  e  (0,T).  This  proves  part  (1)  of  Theorem  (2.9).  Note  that  (B.7) 
is  the  same  as  (2.23)  and,  hence,  we  have  also  established  part  (i)  of 
Corollary  (2.21). 
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UsiriK  Ihr  same  techni  ques,  but  starting  with  (B.tt),  wc>  can  also  prove 
parts  (ii)  of  Theorem  {?.9)  and  Corollary 


^ •  Proof  of  Theorem  (2.13) 

Assumption  (2.  Ill)  means  that  V6  >  0,  i  Tp(6)  >  0  such  that  VT  T^, 

|A^(s)  -  A|  <  6,  Vs  c  (C,1) 

From  (3,3^1  with  |  |*(  • )  |  |^  =  m  and  T  Tp(6) , 

U  S'T.sil  <  exp[ctu(A  ♦  A^(s)  -  A)]  +  (eTm)^  exp(ETm)  (C,2) 

<  exp[cT(u(A)  +  6)]  +  (eTro)*  exp(eTm)  (C.3) 

Now,  choose  as  a  norm  on 

lx|  -  (x'Px)’^^  (C.H) 

where  P  -  P'  >  0  is  the  solution  of  the  Lyapunov  equation, 

A'P  +  PA  +  21  -  0  (C.5) 

nxn 


Under  this  norm,  the  measure  of  a  matrix  M  c  R 


becomes , 


-  ror:  r.  V. 


m(M) 


max  A(M  P 


PM)/? 


(C.6) 


Hence,  (C.3)  becomes, 

l((i(s+T,s)|  <  exp[-cT(1-5)]  +  (cTM)*  exp(ETM)  (C.7) 


By  assumption  (2.1M)  it  is  always  possible  to  select  T„(6)  in  (C.1)  such 
that  6  <  1.  By  inspection  of  (C.7),  these  exists  >  0  such  that  Ve 
e(0,Eo),  |(J)(s+T,s)|  <  1 ,  Vs  e  E^,  which  completes  the  proof  of  part  (i). 

Part  (ii)  can  be  proven  in  an  analagous  manner  starting  with  (B.ii)  and  using 
(C.5)  with  A  replaced  by  -A. 

D.  Proof  of  Theorem  (2.27) 


By  the  variation  of  constants  formula  any  solution  x(s+t),  Vs,t  c  R^, 
of  (2.26)  satisfies. 


S+T 

x(s+-t)  =  4i(s+T,s)x(s)  +  E  j  (ti{s+T,u)  p(u)  du 

s 

where  (J>(s+i,s)  is  the  transition  matrix  for  x  -  e  A(t)x. 

Set 

S+  T 

f(u)  -  J  p(t)dt 

u 

Hence, 

f(s)  -  T  p^(s) 


(D.1) 


(D.2) 


(D.3) 
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where  p^(s)  is  defined  in  (2.28).  Using  (D.2)  and  integrating 
parts  gives, 


x(s+t)  -  4i{s+ T,s)[x(s)  +  eTp^(s)] 
s+  T 

♦  e’  I  (J)(s+T,u)  A(u)  f(u)  du 
s 


Using, 


|<j)(s+-[,u)|  <  exp  (e(s+T-u)]  |a|  1^) 
<  t| |p| 


in  (D.i4)  gives. 


|x(3+t)  -  (Ji(s+T,s)[x(s)  +  eTp^(s)]| 

S+T 

M|A|!.I|PIL  ^  exp[e(3+T-u)|  |A|  |^]du 


<  E 


”I!pIL  (expCetl  |A|  I^)  -  1) 


<  (e  T)* 


exp(ET||A||^)|lp||,,  by  e°-1  <  oe®,  o 


.1)  by 


(D.i*) 


(D.5) 


(D.6) 


q(”.|hlL)||p|L 


(0.7) 


I- 

r 


■ 


which  establishes  (2.?8)-(?.30). 


.‘VV 

•\V>j 


E .  Proof  of  Corollary  (2.31) 


For  any  t  e  ^  integer  k  >  0  such  that  kT  <  t  <  (k+1)T. 


Hence , 


t 

|x(t)|  -  |4i(t,  kT)  x(kT)  +  e  1  ()i(t,T)  p(t)  cIt|  (E.1) 

kT 


<  exp(ETm)( |x(kT) I  +  ex  |1p|1„) 


(E.2) 


where  m  »  11^11^.  Since  (2.32)  holds  by  assumption,  it  follows  that 
|4>(kT,(k-1)T)|  <  exp(-cTB)  where  B  satisfies  (2.35).  Combining  this  with 
(D.7)  for  t«(k+1)T  gives, 


|x(kT)|  <  exp(-eTB)[lx(k-1)T)l  +  eTHp^H^]  +  q(ET  ,m)  ]  [p  |  (E.3) 


c. 


r 


Using  (E.3)  recursively  together  with  (E.2)  and  the  assumption  x(0)  «  0 
gives. 


Ix(kT)|  <  (eT  exp(-eTB)l|p.p||.  ♦  q(ET,m)l  jpj 

•  (1  ♦  exp(-eTB)  +  ...  +  exp(-(k-1)ETB)) 

<  (eT  exp(-ETe)|  |p.j.|  I,  ♦  q(eT,m)|  IpI  |^)/(l-exp(-eTB)) 


•  • » 


*.v 

'.'•’I'-’ 

•  V* 
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.■V  . 


Combining  this  with  (E.D  gives  the  results  (?.33).  (?.3^). 


F .  Proof  of  Lemma  (3.?) 


Using  H(s)  -  d  +  c'(sI-A)  b  we  can  write  (3.1)  as 


e  «  -E(duu'e  +  uc'x),  9(0)  e 


X  “  Ax  +  bu'0,  x(0)  »  0  e  R 


We  now  use  the  Lyapunov  transformation  developed  in  [9]  where 


X  «  Le 


(F.1) 


»  I  ••  % 


m 


(F.2) 


jn  x**.* 


V .  ■:> 


V.'-v 

. .  .  1.  • 


.->■  '.V-  .1 
r»  » • .  •  < 


and  L  satisfies, 


L  -  AL  ♦  bu'  ♦  E(dLuu'  +  Luc'L) 
L(0)  -  0  E  r"’^ 


(F.3) 


We  will  show  subsequently  that  |L(t)|  is  bounded  if  e  is  sufficiently 
small.  In  fact,  if  e  is  not  small  enough  it  is  possible  that  |L(t)|  ->  • 
in  finite  time.  Assume  for  the  moment  that  -3  >  0  such  Ve  e  (0,  e„), 

|L(t)|  is  bounded.  Combining  (F.1)  -  (F.3)  gives. 


e  -  -E  (duu*  ♦  uc'De 


(F.H) 


•  •  •  «  •  •  ^  •  “  «k  •  %  •  •  •  .v^  .  *  *  •  »  ’  « 

«  P  •  »  fc  •  ^  •  k  F  ..  I  ^  ^  ^  k  .  ?  k  >  ^  »  k  ?*  k  F  J 


We  can  decompose  L  as  follows: 


L-Lp+cL,  (F.5) 

where 

L„  -  AL„  +  bu’  (F.6) 

L,  -  AL,  +  d(L„  +  eL,)uu'  +  (L^  +  eLi)uc'(Lo  ♦  eL,)  (F.7) 

Thus,  (F.^)  becomes 

0  =  -e[duu'  +  uc'Lg  +  euc’L,e 

-  -e[u(Hu)'  +  EUc'Lje  (F.8) 

which  is  precisely  the  form  in  (3.3)  with  v  »  c'L,.  It  remains  to  prove 
that  |L(t)|  is  bounded. 

By  the  variation  of  constants  formula,  any  solution  L  of  (F.3) 
satisfies, 

t 

L(t)  -  Lo(t)  ♦el  expr(t-T)A3  GCt,  L(t)]  dT  (F.9) 

O 

where  G(t,L):  x 


is  defined  by, 


C(t,L)  -  dLu(t)u(t)’  ■*  Lu(t)c'L 


(F.10) 


Without  loss  of  generality  choose  c  c  such  that  |c|  »  1.  Then,  V  t  >  0, 


|L|  <  I  ->  |G(t,L)|  <  g()l)|L| 


where 


g(n  -  ||u||jdi|u|L  -  u 


(F.n) 


(F.12) 


■f. 


Since  Re  X(A)  <  0  in  (F.1),  there  exists  positive  constants  K  and  a  such 
that  |exp(tA)|  <  K  exp  (-at),  Vt  c  R^.  For  those  values  of  t  >  0  for  which 
|L(t))  £  I  ve  have  from  (F.11), 

t 

|L{t)|  <  to  ♦  E  k  g(t)  I  exp(-a(t--t))lL(T)|dT  (F.13) 


where  tg  “  ll'-olla,  ^  Bellman-Cronwall  Lemma  [16], 


|L(t)|  <  -_-£KV(t)  ~  exp(-(a  -  eK  g(t))t)]  (F.IH) 


for  all  values  of  t  0  for  which  |L(t)|  <  t,  provided  that  a  -  ek  g(t)  >  0. 
Hence,  |L(t)|  £  I  for  all  t  if 


a/[a  -  i-k  g(t)]  <  !L 


(F.15) 


Combining  (F.lH)  and  (F.15)  gives 


£  <  a  (i-l,)/[K  MhlL{d||u||^  ♦  n] 


(F.16) 


Choosing  I  ~  where 


-  to  [1  *  (1  ^  d||ul|„/tj^^^] 


(F.17) 


maximizes  the  right  hand  side  of  (F.i6).  The  results  in  (3.^)  ~  (3.9) 
follow.  This  proves  Lemma  (3.2). 


0.  Proof  of  Theorem  (3.11) 


Using  Lemma  (3.2)  we  have  that  Ve  e  (0,  e^),  (3.1)  is  equivalent  to 


0  -  -  E  F(t)e  (G.1) 

F(t)  -  u(t)(Hu)(t)'  E  u(t)v(t) 


with  Eg  and  llv||g,  are  given  in  (3.**)”(3.9).  Thus,  we  have  the  local 


average, 


u(t )v  (t ) '  dt 


(G.3) 


-  4  J  u(t)(Hu)(t)’dt  ♦  ^  J 
kT  kT 


(k  +  1  )T 

-  fL.(k)  ♦  4  ^  u(t)[(Hu)(t)-(Hu.  )^(t)]’dt 

‘  kT 

(k>1)T 

+  ^  J  u(t)v(t)'dt  (G.il) 

kT 


where  we  have  used  the  expression  for  R^(k)  in  (3. 18)  and  the  definition  of 
(3.19).  Using  assumption  (3.12)  we  will  show  below  that 


(k+1  )T 

if  u(t)[(Hu)(t)  -  (Hu),(t)]’dtl  <  2(a/b*)| |u| 1*^  (G.5) 

kT 


Hence,  from  (G.i)), 


li[-F^(k)]  <  y[-R^(k)]  *  2(a/b*)|lu|lyT  +  e  ||u|l,||v||.  (G.6) 

<  -  “  *  E  l|u|LI!vIL  {G-7) 


by  (3.15).  Using  Theorem  (2.9)  with  e  <  n/(| |u| |^| |v| |^)  establishes  part 
(1).  Part  (ii)  follows  by  using  (3.17)  to  overbound  u[F,j,(k)]. 

We  now  establish  the  bound  in  (G.5).  For  kT  <  t  <  (k+1)T, 


t-kT  t-rT 

I  h(T)u  (t-T)dT  *  I  1  h(T)u  (t-T)dT 

0  r-0  t-(r+1)T 


by  definition  (3.19).  Thus,  from  (3.12)  we  have 


I  « 


V 

(Hu)(t)  -  1  h(t-T)u(i)di 


k-1  (r+1)T 


I  h(t-T)u  (T)dT  *  I  1  h(t-T 


r«0  rT 


)U^(T)dT 


(by  definition  (3.17)  for  U|^(t)) 


C'''» 


*-.v4 


(Hu)^(t)  -  1  h('t)u  (t-T)dT  *  I  1  h(T)u  (t-T)dT 

t-kT  r-O  t-(r+1)T 


.N*  ‘.V.! 


|(Hu)(t)  -  (Hu).(t)l  <  Mull,  (a/b)  *  I  (g-b(t-(rM  )T)_g-b(t-rt) 

r«0 


Using  this  expression  in  the  left  hand  side  of  (G.5)  gives  the  upper  bound. 


|u(t)l-|(Hu)(t)  -  (Hu),(t)|dt  <  ( |u|  I*(a/b*)(1-e"‘’’”') 


.Cl  .  (e^’^-D  r 


* .  •  .  “  «  ■*. 

*  ■  *  »  •  i«  •  «  •  _  •  •  ^  •  ..  ♦  *1 


.*■  ".***'*>**•'  •*'  •*  -*•  I 


1 


|u||-  1  ^ 


-(k+1 )bT 


1  -  e 


(by  the  geometric  series  formula) 


?  I  |ul  !’(s/bM(l-e 


<  2l|uIP(a/b’) 


which  establishes  (G.5)  and  completes  the  proof. 


H.  Proof  of  Lemma  (5.2) 


The  transition  matrix  for  (5.1)  is  given  by 


S+  T  1 


<ii(s+t,  s) 


(I  +  E  A(t)) 


(H.1) 


S+t-1  S+T-1  * 

I  +  E  I  A(t)  +  E-  I  A(t,)  I  A{t,)  ♦ 

tr*S  tj*S 


+  e^A(s+T-1)  .  .  .  A(s) 


I  +  EtA^(s)  +  R{s,  £■>) 


(H.2) 


by  definitions  (5.3)  and  (5.**).  Hence,  using  the  binomial  series  formulae 
we  have. 


lues.  c,)|  <  (t||A|L)’ . 


i 


V  • 
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